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FINAL REPORT NGR-05-003-UU7 
ABSTRACT 

The mechanism by which atmospheric moisture affects the conductivity and 
dielectric constant of rock specimens is studied in time and frequency domains. 

It is suggested that adsorbed water molecules alter the surface conductivity 
in a manner similar to that observed in semiconductors and insulators. Powdered 
basalts show a low-frequency dispersion produced by the atmospheric moisture 
remaining in the pore system of the sample in a high vacuum; this effect is 

attributed to isolated adsorption centers of the type mentioned above. 

-8 

Simulated lunar permafrost at 100°K and a vacuum of 10 torr together 
with data on lunar samples contaminated with atmospheric moisture and the 
dielectric properties of ice at various temperatures indicate that, if permafrost 
exists in the moon it should present a relaxation peak at approximately 300 Hz; 
for temperatures up to 263°K it may go up to 20 KHz. 

The response of a lunar powder simulator between 100° and 373°K suggests 
that the dielectric properties of the uppermost layer of the lunar regolith 
undergo rapid changes during three tenths of a lunation while remaining 
approximately constant during the remaining seven tenths. 

When a current flows through a rock, interfacial effects arise between 
adjoining mineral grains of different electrical properties; it is concluded 
that in order to have electrical steady state conditions in rock samples it is 
necessary to have volume charge accumulations at interfaces within the sample 
and at the electrode sample interface. A method for measuring heterogeneous 
dielectrics with non-negligible ohmic and dielectric conductivities is proposed 
and experimentally verified. 



2 


Evidence of electrical non-linear behavior in some dry resistive rocks 
is presented, showing the existence of a non-linear region followed by a linear 
one when the voltage, or current density, is increased. The effect is attributed 
to interfacial phenomena occurring at electrode-sample interfaces. A warning 
is made regarding the use of resistivity vs current density plots to decide on 
non-linear behavior. It is suggested that interfacial phenomena and non- 
linear effects may be common features in two-and four-electrode systems. 



INTRODUCTION 


The objective of the .proposed study was to analyze the electrical 
conductivity and dielectric constant of rocks devoid of water, to evaluate 
the changes induced by atmospheric moisture on such properties, to reduce 
or avoid electrode contributions to dielectric permittivity measurements and 
to analyze electrical non-linear behavior in rocks. 

This report is divided in five chapters, each one containing its own 
abstract, conclusions^ references , and figures, and dealing with specific 
aspects of the research items mentioned above. Two chapters consist of re- 
prints of material recently published and the others are either in press or 
being considered for publication. 
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Effects of Atmospheric Moisture on Rock Resistivity 

Roman Alvarez 

Engineering Geoscience, University of California, Berkeley 
Berkeley, California 94720 

The role of water in the electrical resistivity of rocks has long been recognized as an 
important one This study examines the changes in resistivity of rock samples as induced 
by atmospheric moisture An analysis of the term dry rock leads to the conclusion that 
oven drying and evacuation are methods that yield only relative degrees of desiccation To 
obtain absolute values of the moisture present m a sample, one should probably use nuclear 
magnetic resonance and mass spectrometry Experiments were performed on samples of 
hematitic sandstone, pynte, and galena The sandstone underwent a change in resistivity 
of 4 orders of magnitude when it was measured in a vacuum of 5 X 10"’ torr and in air 
of 37% relative humidity Pyrite and galena showed no variations in resistivity when they 
were measured under the same conditions These results, plus others obtained elsewhere, 
indicate that rocks of the resistive type are affected m their electrical properties by 
atmospheric moisture, whereas rocks of the conductive type are not The experimental 
evidence obtained is difficult to reconcile with a model of aqueous electrolytic conduction 
on the sample surface, it is instead suggested that adsorbed water molecules alter the sur- 
face resistivity in a manner similar to that observed in semiconductors and insulators 


Natural state rocks are found m many regions 
to be saturated with aqueous solutions from the 
water table down to a depth of 4 or 5 km 
[Brace, 1971]. At greater depths it is still not 
known how common the pore fluids are m the 
crustal environment According to Wylhe [1971], 
the same mineral assemblages can exist with 
and without a pore fluid , thus the possibility of 
having some aqueous solutions in pores at great 
depths seems indeed to exist 

Above the water table the amount of water 
present at a given time must be a function of 
the local rainfall regime and, to some extent, 
of the atmospheric conditions prevailing at the 
surface The presence of water in the uppermost 
crustal layers of the earth alters in various 
ways the electrical properties of the rocks in 
such regions 

Many studies have been made analyzing from 
various points of view the influence of water 
in rock resistivity [ Hill and Milburn, 1956, 
Mandel et al , 1957 , Keller and Licastro, 1959 , 
Wore, 1964, Parkhomenko, 1967, Scott et al , 
1967, Brace and Orange, 1968], although few 
have dealt directly with atmospheric moisture 
effects Determining the magnitude of the al- 
teration experienced by electrical rock proper- 

Copyright © 1973 by the American Geophysical Union 


ties as a function of moisture content is relevant 
to field measurements as well as to laboratory 
determinations The term moisture shall be 
understood to include any water molecules that 
do not form a part of the mineral constituents 
of a rock sample 

' There seems to be a central problem in trying 
to relate the amount of moisture present in a 
rock sample to the electrical effect it produces, 
the problem consists in precisely determining 
the amount of water in the sample The study 
reported by Wu [1964] in which, by means of 
nuclear magnetic resonance, he was able to 
detect three types of protons associated with 
water molecules corresponding to three different 
types of bonding in clay minerals is a good 
example, given the complexity of these minerals, 
of the nontriviality of the problem Although 
determination of the amount of moisture in 
rocks is mostly done by weighing, it has to be 
recognized that it is only a crude indicator of 
the presence of water in a sample 

Dry Rocks 

When the label dry rock is applied to sam- 
ples used in laboratory measurements, it usually 
has one of three operational meanmgs (1) 
water has not been intentionally added to the 
rock, (2) the rock has been subject to desicca- 
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tion by some heat treatment, or (3) the rock 
has been placed in an evacuated space The 
first two meanings are found most often, the 
second possibly being the most common Few 
instances of the third exist [ Keller and Licastro, 
1959, Howell and Licastro, 1961, Saint-Amant 
and Strangway, 1970], and the vacuums in- 
volved have usually been low (i e , ~10' 3 torr) 

A brief discussion about the water content in 
each of these three cases shall be given 

The sample in case 1 is exposed to atmo- 
spheric moisture, and therefore its surface will 
contain an amount of moisture proportional to £ 
the relative humidity of the air, the interior of £ 
the sample may present a different amount of 
water, depending on its size, porosity, and his- § 
tory g 

In case 2 Worz [1964] has arrived at the a> 

i/l 

conclusion that samples of different thicknesses 2 
dry out differently in an oven, the thicker the c 
sample,' the less the loss of water in a given c 
time at a given temperature Furthermore, this E 
method of desiccation is clearly a relative one, £ 
since the criterion is usually that of observing 
the loss, in weight of a sample, at a given tern- £ 
perature T ,, until no changes in it occur, but, | 
if the temperature is changed to T 2 > T u it “ 
will generally be observed that the previously » 
stabilized weight will be further decreased In M 
Table 1 we show the results of such an experi- 5 
ment in five different rock samples at tempera- " 
tures in the range of those usually reported The ° 
- samples were kept immersed m tap water for 7 ;> 

days, and thereby some degree of saturation is 8 

achieved, whether it was total was irrelevant c 
to the purpose of our experiment The samples 
were allowed to dry (as visually estimated) at " H 
the surface before the first measurement (i e , ^ 

t = 0) was performed The accuracy of the H 
measurements is ±0 0005 gram The oven used 
was of the forced air circulation type With the 
exception of the muscovite schist, all samples 
experienced an extra loss of weight when the 
temperature was raised from 106° to 186°C 
The difficulty in deciding at which temperature 
a particular sample will be considered com- 
pletely dry is obvious Furthermore, if the tem- 
perature is raised sufficiently for a given sam- 
ple, hydration water may begin to be lost from 
its mineral constituents 
Finally, in case 3, moisture near the sample 
surface will be rapidly removed, but moisture 
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Accuracy of the measurements is +0 0005 gram 

♦Although all measurements were completed m less than 2 min, the sample in this one was exposed to room conditions for 10 to 15 min before 
weighing It is likely that the sample absorbed 5 mg of atmospheric moisture in that time 
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in its interior may require long times of pump- 
ing owing to the effect of the pore system of 
the rock, and, if the pores do not interconnect, 
the water may never come out The question 
again is how much moisture is removed in a 
given time and a given vacuum As was pointed 
out above, this manner of desiccating has in- 
volved vacuums of the order of 10~ 3 torr, which 
are usually obtained with mechanical vacuum 
pumps and pumping times of around 10 hours 
In the experiment to be described, evidence is 
shown of water release from the sample oc- 
curring even after 50 hours of pumping at 
around 10" 7 torr Therefore what may seem to 
be a dry rock in a low-vacuum system would 
appear as a strongly contaminated sample in a 
high-vacuum system This phenomenon indi- 
cates that evacuation is also a relative method 
of rock desiccation 

The term dry rock used with the meanings 
described above could be used without compli- 
cations m electrical studies if the effects of small 
amounts of water were in proportion to the 
amount of water present Unfortunately, very 
small amounts of water can have very large 
electncal effects Evidence will be presented 
suggesting that there is a large class of resistive 
rocks that are affected in such properties by 
quantities of water amounting possibly to a few 
monolayers of H 2 0 in the surfaces of the 
sample Consequently, a reconsideration of the 
term dry rock as well as a precise method of 
determining moisture content in rocks seems to 
be needed 

Experimental Method and Results 

To evaluate the effects of atmospheric mois- 
ture on rock resistivity, we used an experimental 
setup consisting of two electrometers, a dc 
regulated voltage supply, a chart recorder, a 
high-vacuum system, and a sample holder The 
electrometers were Keithley models 610B and 
601 having input impedances of 10“ n, they 
were used in two measuring configurations In 
Figure la an electrometer is used to measure 
the total resistance of the sample, it was used 
to monitor the resistance for continuous periods 
of several hours each In Figure 16 the two 
electrometers measure voltage and current 
across the sample The responses of both con- 
figurations were recorded in a Brush Mark 200 
recorder (not shown in the figures), model RF- 



Fig 1. (a) An electrometer R measures the 
total resistance of sample S (6) Two electrome- 
ters, I and V, measure the current and voltage in 
sample S When switch W is in position 1, the plus 
terminal of the voltage source is connected to the 
sample, and dc current and voltage can be deter- 
mined , if the switch is now connected to position 
2, the sample terminals will be open, and a tran- 
sient voltage decay, controlled by the sample 
properties, can be measured If the switch goes 
from position 1 to position 3, the sample termi- 
nals will be short circuited, and the internal cur- 
rent discharge in the sample can be determined 

1783-60 These responses will be further dis- 
cussed The regulated dc voltage source was a 
Princeton Applied Research model TC-602CR 
with a voltage range of 0-60 volts, a no-load 
to full-load regulation of better than 0 0001%, 
and a npple of less than 50-/uv rms at any load , 
The high-vacuum system consisted of a 
Davis and Wilder model 616 system with a 
diffusion pump capable of reaching a vacuum 
of 1 0 X 10' 7 torr m 15 min starting from 
atmospheric pressure, the ultimate vacuum at- 
tainable by the system being in the lower range 
of 10~ s torr The low- and high-vacuum meters 
were those of a Granville-Phillips controller 
unit model 260016 with ranges 10M0' 3 and 
10"‘-10' e torr, respectively The sample holder 
was of the simplest kind and is shown in Fig- 
ure 2 Indium-mercury amalgam electrodes were 
used to ensure an intimate contact with the 
sample, connections to the sample electrodes 
were made through two copper screws mounted 
on acrylic blocks so that high insulation (le, 
>10“ fi) between them and the base plate is 
achieved Electrical connections from atmo- 
spheric pressure to high vacuum were made by 
means of vacuum feedthroughs of the octal pin 
type, the resistance between the measuring pins, 
including as well that of the coaxial cables used, 
was found to be in the 10“-n range, coaxial 
teflon-insulated cable was used for the connec- 
tions inside the high-vacuum chamber 
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Fig. 2. The sample holder used in the measurements consisted simply of copper screws 
mounted on acrylic blocks. The electrical shield has been removed to show the arrangement 
of elements; two samples of different geometries are mounted in place. The copper screws 
make direct contacts with the indium-mercury amalgam electrodes. 


The types of responses obtained by measur- 
ing with the configuration in Figure 16 are 
shown in Figure 3; 1 corresponds to position 
1 of switch W indicated in Figure 16, 2 indi- 
cates switching from 1 to 2, and 3 indicates 
switching from 1 to 3. Let us assume that the 
switch is placed in position 1; the voltage across 
the sample will be the voltage of the source 
(it was kept fixed at +5.000 volts throughout 
the experiments using the configuration in Fig- 
ure 16) less the voltage drop across the current 
meter I, which was kept, at all times, at less 
than 0.005 volt. The correspopding schematic 
response is shown in Figure 3; the response of 
the current to the step voltage applied indicates 
the presence of the Maxwell-Wagner effect in 



TIME 


2 I 3 

Fig. 3. Types of responses obtained for various 
positions of switch W in Figure 16. Number 1 
corresponds to position 1, 2 indicates switching 
from 1 to 2, and 3 indicates switching from 1 to 3. 


the electrode-sample system [ Alvarez , 1972]. 
If the switch is now connected to position 2, 
the sample terminals will be open and the cur- 
rent will go to zero, whereas the voltage will 
decay in a way resembling an exponential func- 
tion. If, instead of switching from 1 to 2, one 
switches from 1 to 3, the voltage will go to zero, 
since the sample terminals will be short cir- 
cuited, and one will be able to measure the 
current flow across the sample; notice the cur- 
rent reversal and decay typical of the Maxwell- 
Wagner effect. The resistance value obtained 
directly with one electrometer was compared 
almost in every instance with the value of the 
ratio obtained with the second configuration. 
They agreed within 5% at all times. The read- 
ings in such cases were made when all transients 
had disappeared. 

A sample of hematitic sandstone (Table 2) 
from the core used in the desiccation experi- 
ment (Table 1) and having a porosity of ap- 
proximately 2.6% was prepared with the in- 
dium-mercury electrodes; the sample had been 
at atmospheric conditions for more than a year. 
Equilibrium with atmospheric moisture through- 
out the interconnected pores of the sample at 
the time of measurement was assumed. The 
sample diameter was 3.22 cm, and its average 
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Rock 

Porosity, 

% 

Resistivity, 
ft m 

Mineral Composition,* 
wt % 


Comments 

Hematitic sandstone 

2 6 

See text 

S3 8 Si0 2 , 26 S Fe 2 0 3 , 11 7 A1 2 0 3 , 

6 8 Ti0 2 

Fine grained 

Pyrite 

0 7 

1 11 x 10- 2 

71 0 FeS 2 , 25 2 CuFeS 2 , 2 8 PbS, 1 

2 Si0 2 

FeS 2 crystal 






size <0 5 ram 

Galena 

0 4 

3 22 x lO" 2 

87 5 PbS, 12 5 Si0 2 


Quartz inclusions 






in massive galena 


Hematitic sandstone from Copper Mine River, Canada Pyrite and galena of unknown origin 
•Determined by X ray analysis 


thickness 8 21 mm , thus the electrodes covered 
an area of 8 13 cm 2 

At the start of the experiment the sample 
resistance at atmospheric pressure was 37 x 
10 s Cl, corresponding to a resistivity of 3 6 X 10* 
Q m The bell jar was located in place, and 
pumping started while the resistance was con- 
tinuously recorded In Figure 4, trace 1 shows 
the sample resistance as a function of time, it 
is clearly seen that the effect of evacuation is 
immediately reflected in a change in resistance 
The initial rate of change of resistivity was 
about 2 1 X 10* Cl m/min Evacuation pro- 
ceeded for 63 hours at pressures from 10‘" to 
10"’ torr At around t = 50 hours a tendency 
of the resistance value to stabilize began to 
appear, to make sure that the readings were 
effectively stabilized, we kept monitoring for 
13 hours more Various rates of resistance in- 
crease were observed while we pumped down 
to the final vacuum attained, toward the end 
of this part of the experiment, variations of 1% 
in resistance occurred in 2-hour periods 


PUMPING STARTED 



30 40 50 60 

TIME (SEC) 


Fig 4 Trace 1 shows the increase in sample 
resistance when the initial evacuation started, the 
resistivity rate of change is about 21 x 10* It m/ 
min After the initial evacuation, air was admitted 
into the high-vacuum chamber, in 23 hours of 
exposure to atmospheric moisture the sample 
reached a value of resistivity of 6.23 x 10* fl m 
Trace 2 shpws the resistance variation as the 
second evacuation' started The resistivity rate of 
change is the same as that of the initial evacua- 
tion 


At t = 63 hours, air of 37% relative humid- 
ity was admitted m the high-vacuum chamber, 
the response of the resistance variation is shown 
in Figure 5 Prior to air admission the resistance 
was stabilized at a value of 2 78 X 10“ Q with a 
corresponding resistivity of 2 75 X 10 s Q m at 
a pressure of 50 X 10' 7 torr, therefore a 
change in resistivity of 4 orders of magnitude 
took place from initial to final conditions Fig- 
ure 5 shows the immediate change in resistance 
when air was admitted, the irregular variation 
around 60 sec corresponds to the lifting of the 
bell jar, so that the sample is directly exposed 
to atmospheric conditions The initial air ad- 
mission is made through a valve, and, once the 
inside pressure is stabilized with the atmo- 
spheric pressure, the bell jar can be lifted The 
transient behavior when the bell jar is lifted is 
systemieally observed We believe that transient 
air currents are forced into the chamber at the 
time of lifting the bell jar and produce the 
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Fig 5. 

TIME (SEC) 

Prior to air admission (1 e , t = 0) the 

resistance was stable at a value of 2 78 X 10" O, 


with a corresponding resistivity of 2 75 X 10*fl m 
and a pressure of 50 X 10~ 7 torr, the sample had 
been in a vacuum for 63 hours At t = 0, air of 
37% relative humidity was admitted, the sample 
resistance started decreasing within 1 sec at a rate 
of 121 x 10 7 n m/min The irregular variations 
around 60 sec correspond to lifting the bell jar to 
allow direct exposure of the sample to atmospheric 
conditions Before the bell jar was lifted, air was 
admitted through a valve to the high-vacuum 
chamber 
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observed response The behavior normalizes im- 
mediately afterward, and the downward trend 
of the resistivity value continues at a rate of 
1 21 X 10 7 Cl m/mm, after 23 hours of exposure 
to atmospheric moisture the sample attained a 
value of 6 23 x 10* Cl m in resistivity The 
• small difference in resistivity between this value 
and the one obtained at t = 0 indicates that 
conditions at this time were very close to those 
at the time the experiment was initiated Exact 
reproduction of the initial resistivity value was 
not attained, probably owing to a 4% decrease 
in the room’s relative humidity 
To test whether the atmospheric moisture 
was, in effect, responsible for the variations 
observed, it was decided to repeat the cycle of 
evacuation and gas admission, except that this 
time dry nitrogen (1 e , 99 99% N 2 ) would be 
admitted instead of air Reevacuation started 
at t = 86 hours, and, after 22 hours of pump- 
ing, the sample reached a value of resistivity 
of 1 03 X 10 8 fimata pressure of 6 4 X 10“ 7 
torr, the response of the rock resistance is 
shown as trace 2 in Figure 4 Notice the same 
rate of resistance increase as in the original 
pumping Further pumping might have taken 
the value of p closer to the one reported pre- 
viously at this pressure range (le, 2 75 X 10 8 
Q m), but it was not considered necessary to 
do so Notice that these values were attained in 
about one third of the time originally needed 
for reaching similar conditions, this reduction 
is possibly due to access paths for the gas/water 
molecules being left open after the first evacua- 
tion, but at present there is no way of verifying 
this assumption The nitrogen was admitted at 
this time until a pressure of I atm was reached 
inside the vacuum chamber The value of the 
resistance was monitored continuously for 
around 2 hours, in which time no significant 
changes of p were observed Finally, a new 
cycle of evacuation and air admission resulted 
in a response similar to the one in Figure 5 
From the observations presented up to this 
point it is clearly established that the reported 
changes in resistivity are due to water effects, 
since, when dry nitrogen substitutes for air, no 
such changes occur Furthermore, we have often 
observed that variations in atmosphenc mois- 
ture (eg, from, 20 to 30% relative humidity) 
correspond to changes in p up to 1 order of 
magnitude 


Our observations have been made on a sample 
of hematitic sandstone, but one can generalize 
the statement made above to a larger group of 
rocks, recalling the results of Howell and 
Licastro [1961] as well as those of Worz [1964] , 
both sets of results have shown the variations 
of electncal properties of rocks as a function 
of moisture by comparing results obtained with 
atmospheric moisture present with those made 
at moderate vacuums (le, ~10“ 3 torr) The 
study of Howell and Licastro [1961] on 23 
common minerals and 71 rock samples concludes 
that below 1 megacycle the presence of moisture 
in a rock plays a dominating role in determining 
its electncal properties, whereas the study of 
Worz [1964] on red sandstone concludes that 
there is a gradual decrease in the frequency 
dependence of c with the decrease in water 
content, whereas for a only a change of the 
absolute value is observed, the frequency de- 
pendence being largely constant Brace et al 
[1965] analyzed the electncal resistivity of 
rocks under pressure and suggested that even 
minute water films may cause the resistivity to 
remain below the dry values by orders of mag- 
nitude 

Consequently, a large number*, of rocks are 
substantially affected by moisture, the magni- 
tude of the change in electrical properties de- 
pending on the specific rock type considered 
There is, however, another group of rocks that 
do not seem to be affected m their resistivity 
by atmospheric moisture We have tried to 
detect changes m the resistivity of samples 
(Table 2) containing amounts of more than 
50% pyrite and others containing more than 
75% galena (presenting p values of 1 11 X 10' 2 
and 3 22 X 10" 2 Cl m, respectively) when they 
are evacuated to ICh 7 torr, and none have been 
observed Brace and Orange [1968] have ob- 
served a similar effect for a Nahant gabbro 
The fact that their resistivities show no change 
when they are measured from atmospheric mois- 
ture to such vacuums indicates that the mecha- 
nism of conduction responsible for their low 
resistivities (rocks are considered of low re- 
sistivity, or sometimes also called 'good con- 
ductors,’ when p < 1 Cl m), whatever it may 
be, has to be considered as dominating over 
atmospheric moisture contributions to conduc- 
tion 
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Transient Responses 

To gam more insight about how atmospheric 
moisture affects conduction in the hematitic 
sandstone sample studied, we shall next present 
the transient responses of current and voltage as 
measured with the arrangement shown in Fig- 
ure 16 In Figure 6a, b, c, and d the current 
and voltage responses of the sample to. the ex- 
citation voltage are shown The time scale in 
Figure 6 b applies to Figure 6a and c as well 
The switching in Figure 66 was made by hand, 
and therefore the length of the pulses is not 
the same, this difference is of no consequence 
to our observations, and it is noted only to 
avoid confusing these dc excitations with a 
square wave excitation often used in rock mea- 
surements Furthermore, the applied voltage 
was always positive 

Figure 6a shows the current response of the 
rock when it is in a vacuum of 5 X 10' 7 torr or 
in an atmosphere of dry N 2 (after evacuation), 


whereas Figure 6c shows the current response 
when the sample was in air of 37% relative 
humidity for 22 hours after the first "evacua- 
tion, this response and that obtained prior to 
evacuation were substantially the same 

Comparing the current responses in Figure 
6a and c, one can determine the change in the 
magnitude of the current flowing through the 
sample These results also suggest that there 
is a modification of the shape of the current 
transient, which is better manifested through 
the voltage decay curves shown in Figure 6 d 
Curve 1 (Figure 6 d) corresponds to the condi- 
tions specified for Figure 6a, and curve 2 to 
those of Figure 6c, they were obtained when 
switch W in Figure 16 was changed from posi- 
tion 1 to position 2 Clearly, the relaxation time 
involved in the high-vacuum or N 2 atmospheres 
is longer than the one obtained at atmospheric 
conditions, and the change in shape of the 
current transients is thereby confirmed Such 
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Fig 6 (a) Current responses when the sample was either in a vacuum of 5 X 10~ 7 torr 

or, after evacuation, m dry nitrogen at a pressure of 1 atm (b) Excitation voltages applied 
to the sample, the time scale shown applies also to (a) and (c) (c) Current responses when 
the sample was exposed for 22 hours, after evacuation, to air of 37% relative humidity The 
magnitude and shape of the current responses have changed with respect to those in (a) 

( d ) Voltage decay curves when the excitation voltage is interrupted, curve 1 corresponds 
to the conditions specified in (a) and reflects the electrical properties of the rock matrix, 
whereas curve 2 corresponds to the conditions in (c) and shows the electrical alterations 
experienced by the rock matnx after exposure to atmospheric moisture 
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behavior is in agreement with a decrease in 
sample resistivity when moisture is present 
Furthermore, the modification of the relaxation 
times implies, through the Kramer-Kronig re- 
lations [Landau and Lifshitz, 1960], that the 
frequency responses are modified as well 

Notice that the voltage decay for both curves 
reaches <~0 9 volt within 10 sec, this remaining 
voltage between the ends of the sample usually 
takes several hours to decay to half its value 
when the sample is, in the vacuum, whereas at 
atmospheric conditions the voltage reaches 
~0 2 volt within 10 min and remains around 
that value for more than 24 hours In fact, we 
have not observed it to go to zero even after 
30 hours of measurements The currents asso- 
ciated with these voltages are usually found to 
be in the range 10 T — 10 -10 amp 

If one tries to observe the transient re- 
sponses of the samples containing pynte and 
galena reported previously, it will be seen that 
the voltage responses when one switches (Fig- 
ure 16) between positions 1 and 2 correspond 
to the shape of the excitation voltage (l e , a 
square pulse), no distortion being observed m 
the time scale of Figure 6 d when either atmo- 
spheric or high-vacuum conditions prevail The 
corresponding currents through the sample are 
square pulses also, thus no Maxwell-Wagner 
effect is observed 

On the basis of these observations, only a 
coarse classification of rocks can be made ac- 
cording to the effects produced on them by 
atmospheric moisture - (1) those whose resistivi- 
ties experience changes of 2 orders of mag- 
nitude or more when alternate measurements 
are made at high vacuums (i e , ~10' 7 torr) 
and at atmospheric conditions with more than 
10% relative humidity and (2) those whose 
resistivities are not affected when they are de- 
termined at the two conditions stated above A 
more precise distinction between these two 
groups, although desirable, must wait for more 
experimental data 

Discussion , 

Once it is established that the resistivity of a 
considerable number of rocks is affected by the 
presence of atmospheric moisture, it is con- 
venient to look into the possible mechanism re- 
sponsible for the resistivity variation 


Conduction through an aqueous electrolyte 
will be considered first, since it is a process of 
conduction that is known to predominate when 
rocks have aqueous electrolytic solutions in 
their pore systems [Brace et al , 1965] In the 
present case, though, we are interested in the 
formation of an aqueous electrolyte on the 
surface of the sample when air is admitted to 
the high-vacuum chamber, since one can assume 
that such an electrolyte is not present when the 
sample is evacuated at 10~ 7 torr 
An aqueous electrolyte is formed when ions 
are dissolved m water. Thus the formation of 
an aqueous electrolyte on the sample surface 
implies the building up of a film of free water 
(le, water not adsorbed to the surface) to 
constitute the aqueous medium in which ions 
from the surface can dissolve and move on the 
action of an external electnc field But it is 
known [ Zisman , 1965] that, if the partial pres- 
sure of a vapor is much below the saturation 
pressure, as is the case for water vapor when 
one admits air into the evacuated chamber, 
the adsorbed monolayer will have a low surface 
coverage Therefore in the first seconds of ex- 
posure to air of 37% relative humidity it is un- 
likely that a water monolayer will form on the 
sample surface, even less likely is the accumu- 
lation of enough water to form an aqueous film 
One would have to assume that times longer 
than a few seconds are involved in the film 
formation, but this assumption would imply 
a time lag m the resistance variation in con- 
tradiction with the experimental observations 
of Figure 5, which show an immediate (i e , 
within 1 sec) response to air admission 
Furthermore, if a continuous electrolytic 
path were to be formed, at time t = t 0 after air 
admission, across the sample surface and con- 
necting the metallic electrodes, the change m 
resistivity should be an abrupt one At this 
time conduction through the highly resistive 
bulk of the sample would be bypassed by at 
least one surface electrolytic path presenting a 
resistivity smaller than that of the rock matnx 
•by several orders of magnitude As more elec- 
trolytic paths were formed, the overall resistiv- 
ity would continue to vary in a gradual manner 
But the rate of resistance variation in the ex- 
periments is never an abrupt one (with the 
exception of the transient when one lifts the 
bell jar) 
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A further indirect indication that conduction 
through an aqueous electrolyte is not taking 
place is given by the results of Fripiat et al 
[1965] They studied the surface conductivity 
of montmonllomtes and silicas containing first 
and second adsorbed water layers and con- 
cluded that the charge carriers involved in the 
surface conduction are not the ions in the mono- 
layer but the protons originating from the dis- 
sociation of water molecules, which is enhanced 
by surface electric fields 

From these considerations it is difficult to 
reconcile the experimental evidence with a 
model of aqueous electrolytic conduction at the 
surface, although other types of electrolytic 
conduction are not excluded [ Scanlon , 1957, 
Snow et al , 1965] Rather, this evidence sug- 
gests a mechanism consisting of an alteration 
of the surface conductivity of the rock sample 
by fractions of a water monolayer adsorbing to 
its surface Such a mechanism would account 
for the immediate variation in resistivity upon 
air admission as well as for the long time (le, 
around 24 hours) necessary to reach equilib- 
rium after evacuation or air admission, these 
periods being necessary for adsorption or de- 
sorption of water in the innermost surfaces of 
the sample, which are connected to the exterior 
through the pore system The surface of the 
sample has to be considered not only as the 
external geometrical surface but also as the 
ensemble of all surfaces in which an air-matrix 
interface occurs 

What kind of evidence is there for assuming 
that fractions of a water monolayer can affect 
the surface conductivity of some rocks ? No 
study in the geophysical literature known to 
this author has shown direct evidence of such 
a mechanism But inorganic dielectrics and 
semiconductors have been shown to be strongly 
affected by various amounts of atmospheric 
moisture At present one has to resort to this 
kind of evidence to support the idea of such 
phenomena occurring in rocks, recalling the 
fact that most rocks in the upper crustal layers 
of the earth can be classified either as dielectrics 
or as semiconductors on the basis of the elec- 
trical properties of the minerals forming the 
rock 

In the rest of this discussion we shall present 
a few of the many available examples of how 
water molecules affect some inorganic siliceous 


matenals We have chosen to present these 
examples to stimulate the idea of similar phe- 
nomena occurring in the upper crustal layers, 
where silica accounts for around 58 wt % of 
the average chemical composition [Ronov and 
Yaroshevsky, 1969] 

A large amount of work has been reported in 
the solid state physics literature about the 
surface properties of semiconductors [ Kingston , 
1956, Many et al , 1965], in which it has been 
recognized that adsorbed atoms or molecules 
may exchange charge with the semiconductor 
and thereby induce changes in its surface con- 
ductivity Adsorbed species can act as electron 
donors or acceptors creating additional allowed 
energy levels for the electrons inside the surface 
of the semiconductor and consequently induce 
changes m the number of available charge 
carriers, these changes take place usually in a 
region a few microns deep 

Buck and McKim [1958], analyzing the 
effects of various atmospheres on the conductiv- 
ity of n type silicon, have shown that it is 
possible to change its surface conduction from 
a strong n type condition to a strong p type 
condition passing through an intrinsic condi- 
tion by a series of chemical treatments and 
exposures to oxygen and water vapor atmo- 
spheres, they conclude that, after any given 
chemical treatment, water vapor makes the sur- 
face more n type and oxygen makes it more p 
type 

Nicolhan et al [1971] report that, when 
water is diffused into a Si0 2 film at tempera- 
tures <200°C and electron currents are made 
to flow through it, a negative buildup occurs 
in the oxide, they attribute this behavior to 
the formation of water-related centers acting 
as electron traps Other types of semiconduc- 
tors have also been shown to be influenced by 
the ambient relative humidity [Kingston, 1955] 

The effects of moisture adsorption on the 
surface resistivity of insulators have long been 
observed, Clark [1962] points out that com- 
mercial inorganic insulators experience changes 
of 5 orders of magnitude in surface conductivity 
upon exposure to atmospheric moisture Fur- 
thermore, he states that the decrease m surface 
resistance' is enhanced by a high degree of sur- 
face porosity, a fact of direct consequence to 
the rock specimens normally handled in labora- 
tory measurements Guyer [1944] presents 
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data showing that the region of maximum vari- 
ation for the surface resistance of fused quartz, 
lime glass, and borosilicate glass occurs between 
30 and 80% relative humidity 

Conclusions 

The problems associated with determining 
the degree of dryness of a rock sample have 
been analyzed, and it has been concluded that 
heating or evacuating a sample produces only 
relative measures of its degree of dryness . It 
appears that, in order to obtain an absolute 
measure of the sample’s water content, one will 
have to resort to methods such as nuclear mag- 
netic resonance and possibly mass spectrom- 
etry Correlations between the absolute amount 
of moisture in a sample and its electrical be- 
havior have to be postponed until the former 
can be properly determined At present, one 
can only study the relative variation of resis- 
tivity between the two different measuring con- 
ditions 

The variations of the' electrical resistivity of 
a hema title sandstone sample measured at 
vacuums of around 10"’ torr and at atmospheric 
conditions, as well as results obtained elsewhere 
on a wide variety of rocks at vacuums of 
around 10‘ 3 torr, indicate that the electrical 
behavior of a large group of rocks is affected 
by atmospheric moisture, whereas our experi- 
ments on two samples of low resistivity indicate 
that there is another group not affected by it 

The possibility of aqueous electrolytic con- 
duction on the surface of the hematitie sand- 
stone sample is difficult to reconcile with the 
experimental evidence An alternative process 
is suggested to account for the rock responses, 
it consists of modifying the surface conductivity 
of the sample by adsorption of quantities of 
H 2 0 possibly as small as fractions of a water 
monolayer Examples of modification of sur- 
face conductivity by adsorbed water molecules 
in semiconductors and dielectrics have been 
presented Because a large number of the rocks 
occurring in the upper crustal layers of the 
earth are either semiconductors or insulators, 
it seems plausible to expect them to exhibit 
similar effects of resistivity alteration by the 
presence of adsorbed water molecules • 

The present study is in no way sufficient to 
prove that fractions of a water monolayer are 
m fact altenng the electrical properties of a 


sample surface, but it seems to give a basis for 
putting forward the hypothesis that surface 
conduction in rocks, as induced by adsorption 
of water molecules, may be an important com- 
ponent of the ensemble of electrical conduction 
mechanisms m rocks 

Such a mechanism would be relevant to the 
study of the conductivity of the upper crustal 
layers in the earth and to the study of electri- 
cal conduction in the moon in the case of exis- 
tence of subsurface water as well as to the effects 
of atmospheric contamination of lunar samples 
on earth Finally, it may account, at least m 
part, for the differences observed in laboratory 
determinations of rock resistivities, which often 
amount to several orders of magnitude among 
the same type of rocks 
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Lunar Permafrost: Dielectric Identification 

Abstract. A simulator of lunar permafrost at 100°K exhibits a dielectric relaxa- 
tion centered at approximately 300 hertz . If permafrost exists in the moon between 
100° and 213°K it should present a relaxation peak at approximately 300 hertz. 
For temperatures up to 263° K it may go up-7o 20 kilohertz. 


The controversy over the existence of 
water in the moon (/) will have an 
“in-depth” test when the electromag- 
netic soundings made during Apollo 
17 have been analyzed (2). The re- 
sults of those experiments will have 
to be interpreted in terms of models 
of the lunar regolith and its dielectric 
properties. 

Permafrost is one of the two possibie 
modes of existence of water in the 
moon, the other being liquid water. 
Models of the internal constitution of 
the moon in which the existence of H/) 
is considered concur m locating the 
permafrost layer above the liquid water 
(2. 4) The depth and extent of the 
permafrost layer is a matter of specula- 
tion. If we assume that such models are 
valid, electromagnetic signals would 
have to travel through permafrost be- 
fore detecting liquid water. 

Information on the dielectric proper- 
ties of lunar samples is available at var- 
ious frequency ranges and temperatures 
(5, 6). However, little is known about 
the role of frozen moisture in the dielec- 


tric behavior of these samples at tem- 
peratures below 273°K. The interpre- 
tation of the lunar electromagnetic 
soundings, on the basis of the data 
presently available, may thus leave ambi- 
guities regarding the presence of the 
permafrost layer. This report provides 
information on the typical response to 
electromagnetic excitations that may be 
expected from a permafrost layer in the 
moon. 

Measurements of the dielectric per- 
mittivity («') and loss tangent (tan 8) of 
a lunar simulator (7) have been carried 
out in the frequency range of 30 to 10 5 
hertz, at temperatures varying from 
100° to 373°K and vacuums of around 
7.0 X 10 -8 torr, in an attempt to 
simulate conditions in the upper layers 
of the lunar regolith. A specially de- 
signed guarded-electrode system was 
used. Complete results of these experi- 
ments will be presented elsewhere (8). 
Here I shall restrict the discussion to 
the responses at 100°K. 

In order to establish a basis of refer- 
ence (9), I measured samples of the 


lunar simulator with packing densities 
of 1.80, 2 00, and 2.20 g/cm 3 contain- 
ing initially atmospheric moisture (at 
40 percent relative humidity). Before 
they were cooled to 100°K the samples 
were maintained at room temperature 
and 5.0 X 10~ 8 torr for periods of 2 to 
3 hours, so that part of their initial 
moisture content was released. In Figs. 
1 and 2 the values of k' and tan 8 for 
the samples at 100°K and 7.0 X 10~ 8 
torr are plotted. Curves 1, 2, and 3 
correspond to the densities cited above. 
Curves 4 (Figs. 1 and 2) were obtained 
for the sample of density 2.20 g/cm 3 
after a cycle of heating (to 373 °K) and 
cooling (to 100°K) took place in the 
high vacuum. The differences with re- 
spect to curves 3 are due to additional 
water losses occurring at the higher 
temperature. This result is in agree- 
ment with observations made of re- 
peated heating and cooling cycles at 
other temperatures (8). Curves 1 through 
4, for tan 8, present broad maxima 
between 100 and 300 hertz, while the 
corresponding k' curves decrease slowly 
with increasing frequency 

The sample that yielded curves 4 was 
used to prepare a new sample of den- 
sity 1.90 g/cm 3 , consisting of the lunar 
simulator (87.7 percent by weight) and 
distilled HoO (12.3 percent). The sample 
had a doughy texture. At 100°K and 
7.0 x 10- 8 torr, the values obtained for 



Fig. 1 (left). Dielectric permittivity (relative) of the lunar simulator at 100°K and 7.0 X 10 -8 torr, plotted against frequency. 
Fig. 2 (right). Loss tangent of the lunar simulator at 100°K and 7.0 x I0 -8 torr, plotted against frequency. Refer to the 
text for an explanation of the curve numbers. 
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Fig 3 Temperature 
plotted against the 
frequency at which 
the peaks of dielec- 
tric relaxation occur 
for ice (12), lunar 
simulator plus wa- 
ter, “dry” lunar 
simulator, and two 
lunar samples (5) 


k' and tan S are shown as curves 5 vapor pressure of ice at 120°K has 

(Figs 1 and 2). If the same data are been estimated as 1 4 X 10 -12 torr (//), 

plotted as k" against «' (*" = k' tan S), therefore, the amount of water frozen 

they correspond to a Cole-Cole distribu- at 100°K and 7.0 X 10~ 8 torr will 

tion (10) given by remain in the sample as long as such 

, , conditions are maintained in the labora- 

= tu' + 1 ^ (1) tory. Curves 3, 4, and 5 (Fig 2) quali- 

tatively establish that the sample losses 
with k 0 ' = 4.90, « x ' = 2 57, r = 6 5 X are proportional to their ice con- 

10 -4 second, and a= 39. The symbol tents A quantitative evaluation of the 

k* denotes the complex (relative) dielec- amount of frozen moisture and its ef- 

tnc permittivity, k„' and k x ' are the feet on *' and k" is presently under way 

values of the real part of k* at zero and (5). 

infinite frequency. The imaginary part The existence of lunar permafrost is 
of k* is / = (— l) 1 ' 4 , <d is 2- times expected within the temperature range 
the frequency, r is the generalized relax- of 100° to 273°K (4) The results pre- 

ation time, and « is a parameter that sented here represent the lower tern- 

can vary between 0 and 1. perature limit It is known that the 

These results show that with increas- relaxation peak for ice moves toward 

ing amounts of moisture the dielectric lower frequencies when the temperature 

behavior of the sample, in the frequency decreases (12) In Fig 3 I have plotted 

range of 30 to 10'’ hetrz and at 100°K, the ice temperature against the fre- 

eventually becomes dominated by ice. quency at which the peak of the relaxa- 

The broad peaks observed in curves tion occurs The ellipses represent data 

1 through 4 (Fig. 2) are interpreted as for ice (12); the triangle corresponds 

relaxations occurring in the frozen to the peak in curve 5 (Fig 2) Extra- 
moisture remaining in the sample. The polatioh of the results for ice is in 

corresponding k' curves are in agree- agreement with my results at 100°K. 

ment with the existence of such relaxa- Lunar samples 10020 and 12022 were 
tions. However, the small k' increments contaminated by atmospheric moisture 
with decreasing frequency preclude its when the dielectric measurements were 
consideration as positive evidence. The performed (5). Consequently, relaxa- 



tion phenomena associated with the 
frozen contaminating water should be 
expected The peaks observed for tan 
S in those measurements have been in- 
corporated in Fig 3 The question mark 
by the lunar samples at 100°K indicates 
the possibility of the peak occurring at 
slightly lower frequencies 

In conclusion the evidence presented 
indicates that if permafrost exists in the 
moon between 100° and 213°K it will 
have a dielectric relaxation peak at ap- 
proximately 300 hertz If its tempera- 
ture is between 213° and 263 °K the 
relaxation peak will occur between 300 
hertz and 20 khz. The frequency at 
which the relaxation maximum occurs 
may serve as a crude thermometer. 

Roman Alvarez 
Engineering Geoscience, 

University of California, 

Berkeley 94720 
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ABSTRACT 

The dielectric response of a powdered basalt simulating lunar fines is 
studied in the temperature range found in the lunar surface while maintained 
at high-vacuum conditions. Two concurrent analyses are made: one regarding 

moisture contamination effects and the other related to the response of the 
lunar regolith. The atmospheric moisture remaining in the pore system of the 
sample in the high vacuum produces a low-frequency dispersion; this is attributed 
to the presence of isolated water adsorption centers in which the local con- 
ductivity values are raised with respect to the conductivity of the basalt. 

It Is suggested that a change in the main electric conduction mechanism 
occurs between 300° and 370°K; close to 300°K conduction seems to be dominated 
by adsorbed moisture, while at 370°K the main conduction process appears to be 
controlled by thermally activated carriers. The dielectric response of the 
lunar regolith during a lunation is schemat ica I fy described. It is found that 
during approximately seven tenths of a lunation the dielectric response of a 5 
to 10 cm surface layer should be fairly constant; in the remaining three tenths 
it should undergo rapid changes. Changes in the dielectric properties of the 
surface layer should be controlled by temperature; in the subjacent layers 
such changes are thought to be controlled by the increasing density of the 
rego I i th . 



INTRODUCTION 


Tho temporaluro on tho surface of Ihe lunar rogolith varies from 
approximately I00°K 1o 400°K during tho lunar synodic period. Tho amplitude 
of the temperature wave Is rapidly attentuated wilh depth ( Roble and Hemingway , 
1971; Hoyt et a 1 , 1971) thus only the first few centimeters experience the 

full change in temperature. In Ihe present study a terrestrial basalt 
simulating lunar powder has been subjected to temperature variations in the 
100° to 373°K range while in a high vacuum. The changes in the dielectric 
properties of the lunar simulator have been determined and are thought to 
represent in a qualitative way the corresponding .variations in the uppermost 
layer (i.e. 5 to 10 cm in depth) of the lunar regolith. 

Although various measurements on the dielectric properties of actual 
lunar samples have been made ( Chung et a 1 , 1970, 1971, and 1972; Katsube 
and Coif ett , 1971; Strangway et al , 1972; Gold et a 1 , 1970 and 1971) few have 
gone to temperatures below 273°K, none of these, however, under high-vacuum 
conditions. Similar studies on terrestrial samples ( Strangway , 1969; Saint- 
Amant and Strangway 1970; Campbell and Ulrichs , 1969) also have not combined 
high vacuum and low temperatures. 

In order to reach conditions similar to those in the lunar surface, a 
terrestrial powder has to lose a considerable amount of the atmospheric moisture 
usually adsorbed to its surfaces ( Alvarez , 1973b); in so doing the sample 
experiences changes in its dielectric properties. An analysis of these changes 
yields an insight into the reverse process, namely that of dielectric variations 
arising from moisture contamination, undergone by lunar samples that have been 
in contact with atmospheric moisture ( Stranqway et a I . , 1972). 

Low-temperature determinations are necessary to generate a complete 
theoretical description of the dielectric behavior of lunar samples, as well 
as to assist in the explanation of electrical phenomena occurring in the lunar 
night. For studies on terrestrial rocks low- temperature measurements offer the 
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advantage of freezing any water that may be present in the sample thus effectively 
reducing conduction effects associated with it. 

EXPERIMENTAL TECHNIQUE 

The experimental set-up consisted of a Davies and Wilder high-vacuum 
system model 616, a cryogenic thermometer with a range of l°K to 300°K from 
American Magnetics, a thermistor calibrated in the 77-400°K temperature range, 
a General Radio capacitance measuri-ng assembly model I6I0-B and a guarded- 
electrode system; all are commercial instruments except the last one. 

A detailed description of the guarded-e lectr.ode system will be given 
elsewhere, here we shall briefly outline its characteristics. A vacuum sealed 
stainless steel box (Figure I) is used as a liquid nitrogen (LN) reservoir 
inside the high vacuum chamber. It is externally fed by means of adequate 
feedthroughs. A heavy copper strap is welded to the base of the LN reservoir; 
as long as there is LN in the reservoir the strap maintains the LN temperature. 

At the top of the copper strap are located the guarded-electrode and the 
guard ring. The former is a copper disk held against the strap by copper screws. 
The latter is a copper ring electrically Insulated from the copper strap by a 
teflon layer. The sample, in the case of solid rock specimens, or the sample 
holder in the case of powders, is placed above the guard and guarded electrodes. 
Heat is directly exchanged between the sample and the guarded-electrode, thereby 
cooling the sample in the high vacuum. The third (upper) electrode clamps 
the sample from above; adjustments in its position are made by means of a 
threaded guide. 

tfeat losses by radiation from the LN reservoir and sample are minimized 
by means of several layers of superinsulation (i.e. aluminized mylar sheets) 

i 

that cover the electrode assembly in the evacuated chamber. In order to reach 
the higher temperatures (i.e. 373°K) pressurized air, or nitrogen, is made to 
flow through the reservoir after being heated to the appropriate temperature. 
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A set of measurements in the frequency range of 30 to 10 Hz, at one 
temperature, involved periods of 6 to 8 hours; of these, 3 to 4 hours were 
necessary to reach thermal quasi-equilibrium conditions, and the remaining 
time was dedicated to the actual dielectric measurements. It must be noted 
that with the present system it is not possible to reach complete thermal 
equilibrium, since heat exchange with the sample occurs only at one of its 
ends, namely the one in contact with the guarded electrode. Measurements were 

made at temperatures of 100°, 298°, and 373°K. The powdered sample was packed 

3 2 

to densities of 1.8, 2.0, and 2.2 g/cm in the working volume (area 17.9 cm 

and thickness 0.5 cm) of the sample holder. This- is made out of a good quality 

dielectric material (NEMA G-10) with suitable mechanical characteristics in 

the temperature range of interest. 

The sample was packed in the sample holder at atmospheric conditions and 

set in place in the electrode assembly. Evacuation started slowly with the 

mechanical pump; in approximately 45 minutes the sample was exposed to the 

high vacuum. At this time the vacuum was usually between 6.0 x 10 ^ and 

3.0 x 10 6 torr depending on the packing density of the sample and the atmospheric 

humidity at which the sample was prepared. It was left to outgas in the high 

vacuum for approximately one hour more; in this time the pressure decreased 
—8 -7 

to 3.0 x 10 - 3.0 x 10 torr. Liquid nitrogen was admitted to the LN 

reservoir and cooling of the sample started. 

The temperature of the guarded-electrode reaches 82°K in 20 minutes, after 

-8 -7 

the first filling of the LN reservoir, in a vacuum of 10 - 10 torr. The 

sample grains directly in contact with this electrode are expected, to attain 
such a temperature in approximately the same time. We observed, however, that 
at places within the sample, not directly in contact with the electrode, heat 
propagation was extremely slow under high vacuum conditions. This is possibly 
due to a heat conduction mechanism in which radiation predominates (Cremers, 
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1972); whether this is the actual case cannot be ascertained from our experiments. 
For calibration purposes several temperature determinations were made inside 
the powder sample. It was found that temperature gradients of approximately 
I 00° K/ cm, in the direction of the symmetry axis of the sample holder, existed 
in the lunar simulator of density 1.80 g/cm^ after 90 minutes of the first 

_8 -7 

filling of the LN reservoir when the vacuum was in the 10 - 10 torr range. 

Thus, the time necessary to reach thermal quasi-equilibrium conditions 
appeared to be unduly long. In order to reduce it, we increased the pressure 
in the vacuum chamber by means of dry nitrogen (i.e. 99.99 % N^), since It 
does not induce changes in the electrical response of "resistive" rock specimens 
( Alva rez , 1973b). It was observed that the thermal conduction of the sample 
improved markedly when the pressure was between 0.1 and 1.0 torr. While 
higher pressures only slightly improved thermal conduction in the sample they 
increased considerably the overall thermal losses of the system. With the 
sample at a pressure of 1.0 torr of ^ the times necessary to reach thermal 
quasi-equilibrium were reduced to the 3 to 4 hour periods mentioned above. 

The thermal gradients, in the direction of the axis of the sample holder, were 
of approximately 20°K/cm when the dielectric responses were determined. 

Having the sample in thermal quasi -equi I ibrium conditions the was re- 

—8 -7 

moved by a re-evacuation, usually to a pressure in the 10 to 10 torr range, 
before the dielectric determinations were made. The high temperature (i.e. 373 e K) 
was attained following the same sequence of evacuation, admission, and re- 
evacuation previously described, while heated air was being injected into the 
LN reservoir (Figure I). After a cooling or heating process the sample was 
left in a atmosphere of around 1.0 torr for 8 to 12 hour periods in which 
room temperature was attained. 
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DIELECTRIC BEHAVIOR 

A description of the material used as simulator of core S/N 2013 from 

Apollo 12 appears in Table I. The dielectric response in the vacuum as a 

function of frequency, temperature, and packing densities of 1.8, 2.0, and 

2.2 g/crrP appear in Figures 2, 3, and 4 respectively. The experimental error 

limits have been estimated, for , as +_ .01 at all frequencies and temperatures; 

and for tan 6 as: .001 (at 30 Hz) or better for values of tan <$ lesser 

than .04; +_ .005 (at 30 Hz) or better for values of tan 4 between .04 and .09; 

and +_ .01 (at 30 Hz) or better for yalues of tan 4 greater than .09. Each 
curve corresponds to one of the temperatures 100°; 298° or 373°K; for a given 

density a sequence of determinations started at I00°K, continuing to 298° and 

373°K, then repeating measurements at some of those temperatures. In order to 

distinguish the order in a sequence we use a subindex (A, B, or C) in the 

temperature value. This information is relevant to the correlation of moisture 

losses with dielectric responses. The subindex shall be dropped when making 

generical reference to a temperature. 

The samples were frozen to I00°K^ as described in the previous section; 
from the rates of outgasing we inferred that a large portion of the original 
moisture content of a sample was lost in the high vacuum prior to freezing. At 
I00°K one expects the remaining moisture in the sample to be condensed in the 
form of ice. The dielectric responses of the samples at this temperature have 
been discussed elsewhere (Alvarez, 1973a). When the sample was allowed to thaw and 
reached 298°K^ some extra water losses occurred due to additional outgasing in 
•the high vacuum. Heating the sample to 373°K^ increased the thermal energy of the 
water molecules remaining in the pore system of the sample, raising the pressure and 
forcing additional water to be lost. Therefore, when the sample cooled off to 
room temperature (i.e., 298°Kg) the total moisture content in the pore system 


was less than at 298°K. . 

A 
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A preliminary Idea of how moisture Is affecting the dielectric response of 
the sample is gained by comparing the dispersions at 298°K^ with those at 298°Kgj 
the latter are markedly reduced (Figures 2, 3, and 4), such a reduction being 
attributed to the lesser amount of moisture. After the first cycle of heating, 
a second one (l.e. healing to 373°Kg and cooling to 298°K^, as in Figure 3) 
produced little further changes in the dielectric response of the sample, corres- 
ponding to comparatively smaller water losses during the second heating. Our 
observations agree with others (e.g. Howell and Li castro , i96J; Strangway et al , 
1972) which have shown that small amounts of moisture affect mainly the response 
of the sample at the "lower frequencies" (the range of these lower frequencies 
depends on the amount of moisture in the sample and may go as high as I MHz). 
Notice the small but measurable effect that moisture release has on the sample 
at I00°K: Figure 4 shows measurements at I00°K^ and at IOO°Kg, the latter being 

obtained after a cycle of heating (to 373°K^) and cooling took place. 

One must bear In mind in analyzing these results that the same amount of 
moisture in the pore system of the sample contributes differently to conduction 
and dielectric response at different temperatures. At I00°K the moisture Is 
frozen whereas at room and higher temperatures it will become water or water 
vapor depending on the actual pressure in the inner pores of the sample. The 
dielectric and conduction properties of ice and water at the frequencies analyzed 
are totally different. 

Interpretation of dielectric relaxations is facilitated by plotting the 
data in curves of k " against where k" = k ' f an <5 .In Figure 5 we 

have plotted such curves for densities of 1.8 and 2.0 g/crrf 5 at temperatures of 
298°K^ and 373°K^. The plots correspond to Cole-Davidson distributions 
( Davidson and Cole , 1951) given by 
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K *= 


< ' 


(I + jOJT) 


3 


(I) 


where <* is the complex (relative) dielectric permittivity, ** and are 

0 co 

the real parts of at zero and infinite frequency, w is 2n times the 
frequency, t is a relaxation time, 3 is a parameter that can assume values 

I /2 

between zero and one, and j = (-1) ; the imaginary part of £ i s denoted k" . 

In Table 2 are summarized the parameters of the Cole-Davldson distributions 
obtained for temperatures of 298°K^, 373°K^, and 373°Kg, and the three densities. 

The data in Table 2 suggest that the relaxation time t is controlled 
mainly by the temperature (in the temperature range of 298° to 323°K) a II hough, 
as will be seen, moisture also has an effect on x at 298°K; they show that 
increasing densities yield larger and values, as expected for a given 
temperature, and suggest that increases in the parameter 8 occur for increasing 
temperatures, although the data for the latter observation do not seem to be 
conclusive. In addition, decreasing amounts of moisture produce decreasing 
values in at 298°K and 373°K; comparison of values at 298°K^ and 298°Kg 
(see Tables 2 and 3), and comparison of k’ values at 373°K. and 373°K_ for a 
density of 2.0 g/crr ? (see Table 2, low-frequency relaxation) substantiate this 
observation. In spite of the change in ie’ between 373°K^ and 373°Kg the relaxation 
time* and the parameters do not experience significant variations, suggesting 
that moisture has little effect on them at 373°K. The 30 Hz frequency-points 
(Figure 5) for 373°K^ and densities 1.8 and 2.0 g/cm^ have been interpreted as 
arising from a residual electrode impedance, based on the data plots at 298°K^ 
which clearly manifest this effect ( Dansas et a 1 , 1967). 

The curves corresponding to the three densities at 298°Kg are plotted as 
< " against k 1 in Figure 6. They could not be fitted to Cole-Davidson distributions; 
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instead two Cole-Cole distributions have been fitted to each set of data 
(Cole and Cole , 1941). The analytical expression for such a distribution is 



+ 



I + (Jut)* a 


( 2 ) 


where the symbols correspond to those of equation (I), recalling the fact that 
the relaxation time has different meanings in a Cole-Cole and in a Cole-Davidson 
distribution and that a is a parameter that assumes values between zero and 
one. Equation (2) and the parameters shown In Table 3 fit the data within the 
error limits ( Ax' and A k?' ) given in the same table. The relaxation times 
between the low and high-frequency relaxations differ by three orders of 
magnitude, strongly suggesting different causative mechanisms. The ^ values 
of the high-frequency relaxation at 298°Kg are essentially the same as those 
obtained for the Cble-Davidson distributions at 298°I< A , as expected from moisture 
effects. 

DISCUSSION 

The above observations establish a general pattern for the dielectric 
behavior of the lunar simulator and the effects introduced by varying amounts 
of moisture. The responses at IOO°K show small dispersions, with maximums in 
tan 5 at around 300 Hz, associated with the frozen moisture in the sample 
( Alvarez , 1973a); at this temperature and in the frequency range studied there 
is no evidence that the basalt powder itself presents relaxation phenomena. 

At 298°Kg (Figure 6) it is observed that two independent relaxations appear; 
the high-frequency relaxation is attributed to the dielectric properties of 
the powdered basalt. The low-frequency relaxation is attributed to interactions 
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between the moisture in the sample and the powdered basalt as follows: It has 

been proposed ( Alvarez , 1973b) that adsorbed water molecules Increase the 
surface conductivity of "resistive” rock specimens by creating additional 
allowed energy levels in the surface of the sample; in the present case we 
propose that the low amounts of moisture involved create a series of isolated 
adsorption centers in which the local conductivity value is raised, by the 
same mechanism, with respect to the surroundings. 

This results in a model equivalent to that of conductive grains (e.g. 
spheres) imbeded in a resistive matrix which gives rise to a low-frequency 
relaxation of the Debye type through a Maxwe I l-Wagner effect ( Koops , 1951; 

A Ivarez , 1973c). The low frequency response In the present experiment! actua I ly 
corresponds to a Cole-Cole distribution of relaxation times instead of a 
single relaxation time (i.e. Debye relaxation); such an effect is attributed to 
a distribution of conductivities in the isolated centers and to the various 
Shapes they must acquire. The plausibility of this explanation Is brought about 
by the statistical nature of the adsorption and desorption processes, coupled 
to the variations in dielectric response induced by different shapes of con- 
ductive inclusions ( Si I lars , 1937). 

When larger amounts of moisture, at the same temperature, are present in the 
pore system one would expect an increase in the density of isolated centers of 
higher conductivity. Accordingly, the magnitude of the Maxwe I l-Wagner effect would 
be increased producing larger values as well as larger energy losses. This 
actually results when the distributions at 298°Kg and at 298°K^ are compared 
( i.e., Figures 5 and 6, and Tables 2 and 3); in the latter the amount of moisture 
Is larger. As it increases in magnitude the low-frequency distribution affects 
higher frequencies creating a smoother transition between the low and the high- 
frequency relaxations. The response of the powder becomes overwhelmed by the 
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moisture effects and the two distributions coalesce into a single one 
approximately described by a Cole-Davidson distribution. The extrapolated k' w 
value remains essentially unchanged during the transition, in agreement with 
expected moisture effects. 

The Cole-Davidson distributions at 373°K^ (Figure 5) present the peculiarity 
of having increased their relaxation times with respect to the distributions 
at 298°K^ (see Table 2). It Is well known that when thermally activated 
mechanisms are present they produce Increasing conductivities and decreasing ' 

t 

relaxation times with increasing temperatures (e.g. Saint-Amant and Strangway , 
1970). According to the increase in relaxation time at 373°K one would have 
to rule out the possibility of such a mechanism operating throughout the 
temperature range of 298°K to 373°K, and yet one would have to account for the 
higher conductivities at 373°K, as obtained from a = * ,, we o where 0 represents 
the ohmic plus dielectric conductivities, to Is 2ir times the frequency, 

k” = tan6 and e = 8.85 x 10 12 f/m. 

o 

One could not attribute the Increase in conductivity to water effects; at 
any rate one would expect a smaller density of isolated centers of higher con- 
ductivity at the higher temperature since water molecules should be excited to 
higher energies, reducing the population of adsorbed molecules. An indication 
that energy losses associated with moisture are not contributing significantly 
at 373°K is obtained comparing the and tans curves at 373°K^ and 373° 

(Figure 3); they suggest that either little moisture remains in the sample or 
that its contribution to losses has diminished, in any event the energy losses 
would not be controlled by moisture at 373°K. 

Thus it seems reasonable to assume that a change in the main conduction 
mechanism takes place between 298°K and 373°K. Up to some temperature value, 
above 298°K and below 373°K, the main conduction mechanism would be controlled 
by moisture and from there on a thermally activated mechanism would predominate 
as described elsewhere (e.g. Saint-Amant and Strangway , 1970; Stranoway et a I . , 
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1972; Chung et a I . , 1972). The last mechanism should cause the overall con- 
ductivity of the sample to be Increased, contrasting with the water mechanism 

that would modify only isolated centers. The thermally activated mechanism 

1 

should mask the contribution to energy losses arising from remaining adsorbed 
water molecules and should account for the increased relaxation time of the 
Cole-Davidson distribution at 373°K. Although the present data are not 
sufficient to prove these statements, the consistency of the results in three 
Independent experiments at three different densities indeed. supports them. 

DIELECTRIC VARIATIONS IN THE REGQLITH 

Some inferences can be made regarding the dielectric variations undergone 

by the lunar regolith as a function of temperature and density. These will be, 

of course, based on the present results and should be reviewed whenever similar 

studies on actual samples of the lunar regolith are available. 

For the purpose of analyzing its dielectric variations the regolith can be 

divided Into two layers: a surface layer of around 5-10 cm depth in which the 

dielectric response shall be controlled by surface temperature variations and 

a second layer, beneath the first one and extending down to the basement, in 

which the dielectric properties will be controlled by density. 

The response of the surface layer can be outlined considering two temperature 

regions: one, between 100° and 300°K, in which <' and tan<$ increase slowly 

with increasing temperature (i.e., corresponding to the present observations in which 

ic t increases up to 6$ and tan6 increases up to .07 from the values at I00°K); 

and the other, between 300° and 373°K, showing rapidly increasing values of and 

tan5 with increasing temperature (i.e., corresponding to increases in k' of 
5 

around \ 0 $ at 10 Hz and up to 65 % at 30 Hz and tanS increasing from .04 to .08, 
with respect to the values at 298°K). The limiting value of 300°K is only tentative, 
we believe that it may go up to 325°K. 

In any event, equal temperature increments will produce different dielectric 
variations in the two temperature regions. When these observations are coupled 



to the temperature variations in the lunar surface during a lunation (Robie and 


Hemingway , 1971), a clearly assymmetric response results for the dielectric 
properties of the surface layer: it will be fairly constant, and approximately 
represented by the response at I00°K, during roughly seven tenths of a lunation, 
while undergoing relatively strong variations in the remaining three tenths 
of the lunation. 

The material below the surface layer (i.e. the second layer) has been found 
to vary rapidly in density with depth ( Mitchel I et al , 1972), while the surface 
temperature variations are strongly attenuated ( Robie and Hemingway , 1 97 1 ) . 

Such a situation results in <’ being mostly dependent on density in the lower 
layer. The relative increments in K ' produced by temperature and density in- 
crements in the present experiments, also establish the preponderance of 
density over temperature as the controlling parameter: Figure 7 shows that 

only below I KHz an Increment of approximately 300°K in temperature and an in- 
crement of .4 g/cm 3 In density (i.e., in the density range of 1.8 to 2.2 g/cm 3 ) 
produce similar variations in ** ; tan $ seems to be rather Independent of 
density in the density range analyzed. 

We found that the data at room temperature (I.e., 298°Kg and 298°K^) do 
not follow Rayleigh's mixing formula ( Campbell and Ulrichs , 1969) in the fre- 

5 

quency range 30 to 10 Hz. The difficulties involved in relating the density 
of the powdered material to its dielectric response, when temperature variations 

5 

are Involved at frequencies below 10 , are obvious in Figure 7. 
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CONCLUSIONS 

A powdered basalt simulating materfal of the lunar regol ith was subjected 

to temperature variations similar to’those found in the lunar surface during 

a lunation. In order to obtain the dielectric response under simulated lunar 

—8 —7 

conditions all measurements were made in vacuums of 2.2 x 10 to 4.0 x 10 
torr which, in addition, allowed for an evaluation of moisture effects in the 
powdered sample. 

The dielectric responses at 298° and 373°K were identified as Cole-Cole 
or Cole-Davidson distributions, and their parameter values determined. Two 
Cole-Cole distributions were found for each set of data at 298°Kg, one corres- 
ponding to a high-frequency relaxation (with x in the 10 ^ sec range) 
and another to a low-frequency relaxation (with t in the 10 ^ sec range). 

The former was attributed to the powdered basalt; it was proposed that the 
latter is caused by a Maxwe 1 1 -Wagner effect which originates in isolated water 
adsorption centers; in such centers adsorbed water molecules would raise the con- 
ductivity with respect to the uncontaminated surroundings. The data at 298° 
and 373°K suggested a change in the main conduction mechanism taking place 
between these two temperatures: near 298°K conduction seems to be controlled by 
moisture, while near 373°K it appears to be controlled by thermally excited 
carriers. 

The lunar regol ith was divided into two layers: a surface layer 5 to 10 cm 
deep in which the dielectric response is thought to be controlled by temperature 
variations, and a subjacent layer in which the variations in k’ are thought 
to be controlled by density; in this layer the tan 6 values would probably be 
slowly varying functions of density. An assymmetry in the dielectric properties 
of the surface layer during a lunation was suggested: during approximately 

seven tenths of a lunation this layer would show a slowly varying behavior in 
Its dielectric properties, undergoing rapid variations in the remaining three 

tenths of the lunation. 
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TABLE I. Source and Sizes of Basalt Components of Lunar Simulator* 
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TABLE 2. Parameters of the Cole-Davldson Distributions at 298°K^, 373°K^, and 373°Kg 
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mainly in the frequency range of 100 to 1000 Hz. 



TABLE- 3,- Para.mete.rs. of the two Cole-Cole Relaxation Mechanisms at 298°K 
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LIST OF FIGURE LEGENDS 


Fig. I. Schematic view of the guarded-electrode system. The whole assembly 
is maintained in the high vacuum; low and high temperatures are ob- 
tained when liquid nitrogen fills, or heated air flows through, the 
reservoir. The sample exchanges heat with the reservoir via the 
copper bar and the guarded electrode. 


Fig. 2. Ca) Dielectric permittivity against frequency and (b) loss tangent 

1 5 

against frequency for the sample of packing density 1.8 g/cm . 

Temperatures followed the sequence I00°K^, 298°K^, 373°K^, and 

—8 

298°Kg. The pressures ranged from 2.2 x 10 torr at I00°K^ to 

1.2 x I0~ 7 torr at 373°K.. 

A 

Fig. 3. (a) Dielectric permittivity against frequency and (b) loss tangent 

against frequency for the sample of packing density 2.0 g/cm^. 

Temperatures followed the sequence I00°I< A , 298°I< A , 373°I< A , 298°« B , 

—8 

, and 298°K^. The pressures ranged from 3.0 x 10 torr at 
I00 °K a to 2.5 x I O' 7 torr at 373°^. 


373<>K B 


Fig. 4. (a) Dielectric permittivity against frequency and (b) loss tangent 

against frequency for the sample of packing density 2.2 g/cm 3 

Temperatures fol lowed the sequence I00 0 « A , 298°K A , 373°K A , 298°Kg, and 

-8 

IOO°Kg. The pressures ranged from 2.5 x 10 torr at lOO 0 !^ to 
4.0 x I0" 7 torr at 373 0 K A . 

Fig. 5. k?' against plots for the distributions at temperatures of 298 0 « A 

and 373°« a , and packing densities of 1.8 and 2.0 g/cm . The dis- 
tributions are of the Cole-Davidson type and the sets of parameters 
describing them appear in Table 2. The numbers by the data points in 
the distribution of 373 0 « A and of density 1.8 g/cn? represent fre- 
quencies in kilohertz; to avoid overcrowding only a few frequencies are 


indicated. 



l+o 


Fig. 6. k " against k’ plots for the distributions at a temperature of 298°Kg 
and densities of 1.8, 2.0, and 2.2 g/cm" 1 . Each one is described by 
two Cole-Cole distributions represented by dotted curves; the corresponding 
sets of parameters appear in Table 3. The numbers by the data points 
of the distribution of density 2.0 g/cm^ represent frequencies in 
kilohertz. 

Fig, 7. Dielectric permittivity against density, with temperature and frequency 
as parameters. Increments in density, from 1.8 to 2.2 g/cm^ , pro- 
duce larger increments in K' than those produced by temperature variations 
from 100° to 373°K. Only for frequencies below 1 KHz the temperature 
effects become comparable to density effects in the ranges analyzed. 
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Complex Dielectric Perm iff i vity in Rocks: 

A method for its measurement and analysis 
Rorr^n Alvarez 

Engineering Geoscience 
University of California 
Berkeley, Calif. 94720 

ABSTRACT 

The interfacial effects arising in electrical measuremenls 
of rocks, when a two-electrode system is used, are analyzed via 
the Maxwe I l-Wagner effect. The similarity in electrical behavior 
between rock samples and heterogeneous dielectrics with non- 
negligible ohmic conductivities leads to a general analysis of 
the Maxwe 1 1 -Wagner effect, from which it is concluded lhal in order 
to have electrical steady state conditions in rock samples it is, 
in general, necessary to have volume charge accumulations at inter- 
faces within the sample and at the electrode-sample interface. 

The expression for the charge accumulation at the inlerface as a 
function of frequency is obtained as well as expressions for the 
effective dielectric permittivity in a composite material; it is 
shown tljat the relaxation time associated with the charge accumulation 
and that of the effective dielectric permittivity are the same. 

Given a sample consisting of two different materials it is 
shown that whenever both components have non-negl igible ohmic con- 
ductivities, the imaginary component of the effective dielectric 
permittivity tends to infinity when the frequency tends to zero; 
such behavior is not observed when one of the two conductivities 
is negligible. This fact constitutes the basis of the meihod of 



measurement for heterogeneous dielectrics witli non-negl igibie 
ohmic conductivities proposed herein; it consists essentially of 
utilizing a dielectric with negligible ohmic and dielectric con- 
ductivities as one of the two components of the sample, which 
isolates fhe o^her one from the electrodes. 

Theoretical examples of the behavior of such a composite 
material for various cases are analyzed and three experimental 
examples are presented. The method presents the advantage of 
permitting application of theory developed for dielectrics in 
genera I . 



INTRODUCTION 



Natural rocks are heterogeneous materials of considerable complexity 
from an electrical point of view. A given sample can be thought of as a 
multi-crystalline matrix with a system of pores among which a given degree 
of interconnection exist; these pores may contain some free water and, if 
ions are present in it, an electrolyte; in addition adsorbed water and gases 
are found at the matrix-air interfaces. Each component contributes to the 
electrical properties of the rock. The conductivity O' and the dielectric 


c | c - 

permittivity c are the usually determined properties. 

We shall deal with the problem of conduction effects at the electrode- 

2 

sample interface when a rock sample is placed in a two-electrode system. 

These effects are most clearly manifested in the dielectric permittivity 

measurements of rock samples at frequencies below I KHz, although this figure 

varies according to rock type and conditions of measurement. In any event, 

when the frequency tends to zero hertz, the values of dielectric permittivity 

increase according to an inverse power of frequency; the power value is usually 

found between I and 2, varying with author. The relative dielectric permittivity 

( K = , where £ 0 is the vacuum permittivity) reaches values of around 10^ at 

4 ix 

frequencies of less than 100 Hz. At frequencies above 10 Hz the values of rv 
are found to be "normal", where normal means the values to be expected from the 
contributions of dipolar, ionic and electronic polarizabilities in a given 
material. In the geophysical literature the values of |0 6 have often 

been called "abnormal", the abnormality stemming from the impossibility of ex- 
plaining such values in terms of the mechanisms of polarization cited above. 


MKS units are used in this work. 


2 


For a summary of other systems and their problems, see Fuller and Ward, 1970. 



In homogeneous dielectric substances the value of }{ ai zero frequency seldom 
2 

goes above 10 . 

Rock measurements are not the only insfances ofK-values higher than 10^ 
at frequencies of less than I KHz; they have also been observed in studying 
heterogeneous dielectrics and some types of semiconductors as well as in 
electrolytes. 

Attempts have been made at evaluating or supressing the high H -values 
in rock samples at low frequencies. Some researchers have used mica sheets 
placed between the sample and electrodes (Tarkhov, 1948; Valeev and Parkhomenko, 
1965) thereby achieving considerable reductions in the «. -values. Others 
(Scott, et al, 1967) have rejected this scheme in view of experimental problems 
and the lack of an adequate theoretical model to analyze the results. 

The departures of K.from normal values become accentuated when measuring 
rocks with a considerable amount of moisture. Parkhomenko (1967) analyzes 
this problem and concludes that minor variations in water content may give rise 
to extreme variations in dielectric permittivity. 

For "dry" rocks, and direct electrode-samp le contact, the values of fC 
have been found within reasonable limits for frequencies ^ 100 Hz (Keller, 1966, 
p. 569; Parkhomenko, 1967, pp. 38-41). However, a tendency for 1C to increase 
has been observed at frequencies lower than 100 Hz even in rocks that have been 

oven and vacuum dried, (Keller and Licastro, 1959, p. 263). Furthermore, in- 
creases in K. are systematically accompanied by higher losses (i.e. higher 
energy dissipations). 

The above discussion is based primarily on measurements in rock samples. 
Summarizing those results one can establish that (i) the influence of the direct 

contacts between sample and metallic electrodes gives rise to high apparent values 

of dielectric permittivity K, , (ii) such values can be reduced to normal values 

by avoiding the direct contact between sample and electrodes and (iii) a direct 

relation seems to exist between the ohmic conductivity of the sample and the high 
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values of . 

The similarity between the electrical behavior of rocks and that of 

heterogeneous dielectrics with non-negi i g ib le ohmic conductivities is illustrated 

by the works of Blechschmidt (1938), Koops (1951) and Krotzsch (1964), which also 

3 

show values of tC** ® for frequencies below IKHz accompanied by high energy losses. 
Dieleclric liquids in contact with metallic electrodes yield similar results 
(Johnson and Cole, 1951; Dansas et al, 1967). Surprisingly the analyses made 
for these materials have not been fully applied to conduction phenomena in rocks. 

As pointed out by Johnson and Cole (1951) the combined effects of dc (ohmic) 
conductivity and electrode polarization give rise to, serious errors in the observed 
apparent dielectric permittivity and energy tosses. However, suppressing electrode 
polarization does not suppress the effects arising from the heterogeneity of the 
sample , which will give rise to values of dielectric permittivity and energy 
losses higher than those obtained from homogeneous samples of equivalent compositions. 

THE MAXWEIL-WAGNER EFFECT 

Discussions on the different causes of polarization in dielectrics most 
of the time include descriptions of dipolar, ionic and electronic mechanisms 
contributing to the total polarizability in given frequency ranges (Kittel, 1968; 

Hill et al 1969). A fourth type of contribution arises in heterogeneous materials 
due to accumulation of charges at the interfaces but its description is usually 
omitted possibly because, as Kittel (1968) points out 

"This is of little fundamental interest, but it is of considerable 
practical interest because commercial insulating materials are usually 
heterogeneous . " 

Given the inherently heterogeneous character of rocks it should be recognized 

* 

that, at least from a geophysical point of view, interfacial polarization may play 
an important role and deserves careful consideration in the context of electrical 


conduction in rocks. 



The Maxwel l-Wagner effect relates to the charge build-up at the interface 
between two media differing in conductivity and/or dielectric permittivity when 
a current flows across the interface. This phenomenon was first analyzed by 
Maxwell (1873), sec. 328; later on Wagner (1913) studied the effect of spherical 
particles of given conductivity and dielectric permittivity homogeneously dis- 
tributed in a material of differing dielectric permittivity and zero conductivity. 
He found that the relaxation spectrum of such a material had a single relaxation 
time, presenting what is now known as a Debye relaxation in dielectrics (Debye, 
1929). 

Various treatments of this effect can be found in Hill et al (1969), p. 58 
and p. 282, Von Hippel (1966), p. 228, Zheludev (1971), p. 474, and Arnoult et al 
(1965); they are equivalent inasmuch as they consider the two dielectrics having 
no dispersion and, thus, all the losses are due to the presence of the ohmic 
conductivities. It may be recognized that these conditions represent a 
particular case of a more general treatment, namely that of dispersive dielectrics 
with different ohmic conductivities. 

We have approached the theoretical treatment of the Maxwe I I -Wagner effect 
in a way that differs from previous treatments. We consider it more general 
for it does not impose restrictive conditions in the dielectric permittivities 
of the media involved (i.e. we do not require Kj and constants with frequency) 
neither is it confined to the case of sharply defined interfaces between the 
media. 

Consider equations 

\7-D = 5> J=e£ $ = c£ 

where the symbols correspond to the well known variables of electromagnetic theory; 
they are considered averaged over elements of volume containing several molecules 



of a given material. J). J , and 0 are frequency-domain variables. (T* , 

•*** Km- 

^ , and ^ are frequency dependent. The dielectric permittivity ^ and the 
ohmic conductivity will, in general, be functions of position. 

With the continuity equation in the frequency domain 


V‘2 =y-f 


( 2 ) 


and relations (I) one can write 


+ ^V-6 =? <3, 






Defining 


=j"f 

^ ^ ‘ it 

^ € 'J OJ 


the last expression becomes 

/?•£ 4 £-lfe*=o 


and from Gauss' theorem the condition 
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is satisfied across an interface, being the normal component of 

the interface. 

The treatment of the classical Maxwel I -Wagner effect requires and (N 
constants with frequency. Let these conditions hold in the following discussion. 


e 


to 


»<***< 


Further ahead we shall come back to the case of dispersive dielectrics. From 
equations (2), (3), and (4) 

O' a/,.a . W m)l 

( 7 ) 


n ;^= - £ fvCvfr) • L (aj) l 

J J C •» C. L / J 


and the steady state condition V'2 = 0 (this implies either = 0 for 
a/ 4 0 or Cl) = 0. The case ^ = 0 is trivial; therefore we shall consider = 0) 


p(o}= Q* 



one obtains 

(8) 

where 0(0) and £E(0) indicate the values of P(C^) and £. (OJ) at OJ = 0. 

J J \A~ 

It can be seen that if £ and/or vary as a function of position in a 

given volume, the value of V (%) will be non-zero; consequently charge 
accumulations will occur as defined by in equation (8). If & and (4 are 
uniform throughout such a volume, then ^7=0. The term is the one that 

gives rise to the Maxwel I -Wagner effect. 

Therefore if a material is homogeneous it has no volume charge concentrations 
in the steady state, whereas if the material is inhomogeneous ' — ' 
steady state conditions cannot be realized without the presence of a volume 
charge density as defined by equation (8). Furthermore the volume charge 

0(x, y, z) is a function of the electric field (x, y, z) at that position. 

J »***> 

Equation (8) stresses the fact that in order to have a Maxwel I -Wagner effect 
it is sufficient to have either or O' varying in a region, a fact often obscured 
in the classical presentations by the assumption that the two media differ in both 
£ and (K Notice, however, that no Maxwe I I -Wagner effect will appear if the 
ratio has the same value in both media, even if the individual values of ^ 
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and (?- are different for each medium. 

Another aspect of the Maxwel I -Wagner effect is brought about by the vector 
relation between the gradient of the properties and the electric field: only 

the normal component of the electric field to the interface will contribute to 
space-charge build-up. The interface we refer to is a surface, say z = const, in 
a cartesian coordinate system, that has the same £ and values at any point 

on that surface; a volume in which 6: and/or 0^ vary would be made out of a 
succession of such surfaces. 

Time-domain excitation 

Equation (7) can be written in time domain by means of a Fourier transformat ion: 


W-z 


v(%) -E(t| 


(9) 


where the dependence of p on time has been explicitly indicated, and £ are 
constants in time and E(0 ^ the time dependent electric field. 

By requiring now that T7.J = 0, a time independent relation is obtained 

£ = ^ V (%) ■ E « 


( 10 ) 


Equation (10) holds for dc steady state flow of current in a material in which 
J/q)^o . In particular, i f a dc excitation is applied to a sample consisting 
of two slabs of area A and thicknesses dj and with electrical properties 

Q- 


and ^2' ^2 respectively, a charge concentration will build up at the 


interface between media I and 2. When steady state conditions are reached the charge 
concentration is given by equation (10). 

The charge build up can also be seen through Zheludev's (1971) derivation for 
the current density and electric fields in such a material: 

From the continuity of conduction plus displacement currents at the interface 




+- £ - ■ *- 2 — 


(II) 



6o 


and the condition 

E,(W, 4 E 2 Md 2 = V 

where V is the externally applied dc voltage, he obtains 

r/i\ 2V ^4. 

and ' ~ (Mt A4) (M 2 A4) 

T/i\_ (L^ 3 , WlL 

wj+h (6, da 4 - ^4)^0^^2 + ^z. 4) (^ciz'^df 

7 _ ^dz^^di 

' “MT^AdT 

where vj is given by either side of equation (ll).. Notice that 

T. J = 333— We sWq Afe. 


( 12 ) 

(13) 


(14) 


The response of the current to a voltage V applied between times t = 0 
and t = t j ; is shown in Figure I. At t = 0 the current J acquires its maximum 

value. The first term on equation (14) tends to zero as t increases while the 
second remains constant. The first term is known as the absorption current and 
it flows as long as the accumulation of free charge continues at the interface. 

The second term corresponds to a residual current (i.e. the current flowing 

in the steady state). When the external voltage is supressed at t = tj the charge 
stored at the interface is retrieved in the form of a current described by the 
first term in equation (14), but of opposite polarity. 

This type of response we have systematically observed in resistive rock samples 
(i.e. granite, limestone, etc.) when metallic electrodes are in contact with the 
sample (Alvarez, 1971). Low resistivity samples (i.e. pyrite, galena, etc.) in 
contact with metallic electrodes do not show such a behavior. Similar examples can 
be found in Parkhomenko (1971). 



6i 


The previous development corresponds to a "sharp" Interface between media 
I and 2 (i.e. the electrical properties change abruptly at the interface). In this 
case the volume charge distribution ( ^ ) becomes a surface charge distribution 
( X ) since only at the interfacial plane VO) * ° 

However, an equivalence can be established between the charge stored at a 
sharp interface and the charge stored in a volume (Alvarez, 1972, p. 21). Consider 
two media of ,6^ and » ^2 separated by a volume if , of area A normal to 
the applied field, in which and 0 are arbitrary functions of position . Then, 
under steady state conditions, the amount of charge Q stored in the volume is 

the same as the charge stored at a sharp interface of equal area A between the two 
media, provided the current density Jo is the same in both cases. The surface charge 
density is given by 

I = X - 1 - |) 


With equations (12) and (14) (steady state) substituted in equation (15) the 
surface charge density is obtained in terms of the fields E| and E 2 (steady state 

values) in materials I and 2, as 
Frequency-domain excitation 





We turn now to the analysis of the response of a material presenting spatial 
variations in £ and/or \j~ when the applied voltage is of the form 

V-TLe* 


We are interested, as a first step, in determining the behavior of the charge con- 
centration as a function of frequency at the interface between two materials in which 
£ and O' are frequency independent. 



( 17 ) 
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where ^ (°) is the charge concentration at U-/ = 0 given by equation (8). This is 
the maximum value Q can attain at the interface for a given and a given contrast 
between the values in both media. For OJ y O <p ^ J furthermore if 

There is a phase lag between the applied field and the charge concentration 


given by 


<jb = - To 




corresponding to a loss tangent 


7o» S = 




From these considerations it can be seen that the charge concentration behaves 

in the same fashion, except for some constants, as does a dielectric presenting a 

\ 

Debye relaxation, i, e. such a dielectric has zero ohmic conductivity and is represented 


* 

^ 4- 


I 4- 


where ^ * is the complex dielectric permittivity and a and b are constants. This 

I 

comparison hints immediately at the origin of the Debye- 1 ike behavior of a substance 

j 

presenting the Maxwel l-Wagner effect. 

Equation (17) cannot be directly applied to the case of a sharp interface to 
yield the surface charge density as a function of frequency. The difficulty 

arises from the values to be assigned to ^ and O'* at the interface (Alvarez, 1972, 
p. 27). A possible way to obtain for the two- layer material we have been 

discussing is as follows: 

From equations (II) and (12) in the frequency domain 


c W= (&*>**) V 

1 + 4- j£0 * ^ 1 ^ 2 ) 
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The surface charge density at the sharp interface is given by 
substituting equation (19) in equation (20) 

z’m- z(0) 


( 20 ) 


1 +ju)^ 


( 21 ) 


where 


% 


1 (O) « TTTTa ^ ^ l . ^ 

0| CP z d, tfjd. + tfA 


The symbol sba ^ be used * n analogy with the’ complex dielectric permittivity 

^ * to denote the complex character of the function. The real and imaginary parts 
of will be denoted and respectively. * P lo+s as a circle in 

the complex plane. 

Complex dielectric permittivity: definitions. - The expression 

- * ^ 7 ' - ft 

€ £ " J 6 

usually represents the complex dielectric permittivity of a homogeneous material. The 

-/ 

real part ( e ) corresponds to the dielectric "constant" and the imaginary part 
( £" ) to the energy losses in the material. 

Energy losses arise from two different mechanisms: (I) rotation of dipoles 

or bound charges and (2) motion of free charge carriers. Thus, C contains two 


parts 


// 


£ = £ 


// 

Diet- / 


ft 


OHMtd 

CO 


// 


It is common practice to associate the dielectric losses £ to a "dielectric 

DIBL 


conductivity" 


= <o £ 


II 


$ 

4 • 


£•// ft|EL + 


'blCL ’* ,CL 0L) 

But, unfortunately, lumping and in one single conductivity 

has led to the widespread but erroneous practice of applying it in relations as 

G'-coe.' or Q--uje'1oy)6 
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where O’ , then, is implicitly considered to be °r ^(££_ only. The 

experimentally determined parameters are usually <£ ' and fan(f. 

Under some circumstances the error will be small, that is when 
(or viceversa when it is assumed that there is only ^Diet i n the sample) . But, 
for resistive rock samples C)omhic ma V comparable in some frequency range 
to OoieL and the error will be considerable for that case. 

Warnings about such a way of proceeding have been previously made (Campbell 
and Ulrichs, 1969), but for some reason they have not received adequate attention 
Furthermore, there are standard procedures for evaluation of both contributions 
when simultaneously present in one sample (Brockman and White, 1971). 


DieL 


The method to be introduced allows for the determination of both 


£ 


and (j^UMIC 


DIEL 


The following convention will be adopted: 

6 ** £ 'j ^ 


(M 


where 


£ - e'- i e" 


J ^OIEL 

There Is no need to maintain the subindexes DIEL and OHMIC; it will be under- 
stood that ^ will refer 6nly to and ^ +o * Thus 

e J c J co J OJ (22) 

Finally, whenever the material is not known to be homogeneous for certain, 

misleading values of dielectric permittivity may be found experimentally (Hill et al, 

1969) and the ratio C/C should be considered as a relative effective (or apparent) 

o 

dielectric permittivity K r • where C is the experimentally determined value of 
the capacity of the sample and C q is its geometrical capacitance. Furthermore 

% = % 

where is the free space permittivity and 1S e ff ec ‘*’' ve (° r apparent) 

dielectric permittivity of the sample; 6^ may be complex and in such a case it 
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will be written as &Gcl • Since we know that the samples we are dealing with are 

* 

not homogeneous wi I I be used throughout the following developments. 

* 

General expressions for £ eff Equations (I) through (6) were established for <r 

and frequency dependent. In the subsequent developments we analyzed various 
cases in which £ and were restricted to be frequency independent. We shall 
now consider £ to be frequency dependent while maintaining frequency independent. 
These conditions correspond to dispersive dielectrics having constant ohmic con- 
ductivities. The case in which the ohmic conductivities are frequency dependent 
shall not be treated here, although the same forma I i sm .can be applied in such a case. 


Consider the two-layer material described previously. From equation (6) 

* C C n C __ 'M 

Jh C-pcr - 


€. i £ i 


i^i ajLW " cL 

where \f is the voltage across the sample and ^ is its total thickness. The 
interface is considered normal to the field; thus subindex Vj for the normal com- 
ponent of the fields has been dropped. 

Substituting equation (19) in the above expression one obtains 

((ft 4 4 0^,) 4 ja)(^ ( d 2 -f e 2 d) 

which can be finally written as 

-* _ (d*d 2 ) £> (Lz 

1fd^e*A 


COL 1 


(23) 


(24) 


or as 


* 




(25) 


Equations (24) and (25) are the general expressions for the effective dielectric 
permittivity of a composite material made out of two homogeneous slabs of area A and 
thicknesses d ( and d and characterized electrically by tr j and tr 2 * Equation 
(25) gives the explicit dependence of on the charge concentration at the 
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interface. Furthermore, the relaxation time obtained for 


Z 


* and the relaxation 


time associated with the ctr °f a composite material with non-dispersi ve di- 

TT 

electrics (Von Hippel, 1966; Zheludev, 1971), are the same. 

■4( 

The experimentally determined quantities are the components of £ grr , 
namely tZ and c ^ . The two limiting cases U)~50 and , for 

the components of , are of interest: 

From equation (25) 

€ ' — -> i!L 

e, d 2 + CTVd, 


A)— = > o 


co — > O 


(26) 


£ 


// 


OJ — =>«> 


£r,£ z 

€.\ d z +€zdi 


and £ 


if 


^ O 


CO 


<x> 
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where we have written 6* and ^ instead of ^ and because ^ and will be 


zero when and when <x> ~^ao . 


.// 


From the limiting case (0-^C> it is seen that the absolute value of (f ^.— >~°° ; 
we shall label this a "divergent behavior"; thus, there will be cases in which it 
will become very difficult to determine experimentally, with any accuracy, the two 
components of c - ^ in this limit since one of them will be much larger than the 
other. It is worth noticing that the divergence occurs only for one component and 
only in one limiting case; in fact, the problems arise from the two ohmic conductivities 
( \ and $2 be ing non-negl igi ble . , If one of them (i.e. (j^ in equation (25) is 
negligible, such divergence is eliminated and the two components of C7 ^ become 
finite in the two limiting cases expressed by relations (26) and (27). 

METHOD FOR DETERMINATION OF 6: AND IN ROCK SPECIMENS 
The scheme for measuring electrical rock properties we shall introduce consists 
of placing a "good" (or lossless) dielectric material (i.e. one in which 
and are negligible) between the electrodes and the sample to be 

analyzed. The former is, of course, considered as homogeneous, whereas 
the latter is known or assumed to be heterogeneous... 


The composite sample, in 
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any event, will be heterogeneous. We want to apply the results obtained fon 

two homogeneous materials to a system composed of one homogeneous material and 

one heterogeneous material. This way of proceeding is justified inasmuch as 
the heterogeneous material can be represented by an ef fecti ve dielectric 

permittivity which may itself be constituted by two or more different electrical 

components . 

Therefore we shall be dealing with two Maxwel I -Wagner effects of different 
origi n ; one, deliberately established and controlled bv the interfaces between 
the good dielectric and the sample and, the other, occurring within the sample 
and controlled by its constituents and their particular distribution (i.e. 
geometry). The first one we can alter by varying the thickness of the good 
dielectric as will be shown later, the second we cannot modify without changing 

the properties of the sample. 

c* 

Once C e y has been obtained experimentally as a function of frequency it 
will be possible to model it theoretically; this can be done by assuming a 

c* 

specific functional form for and fitting the model parameters to the data. 

The final properties of the sample will be given as a set of parameters, 

one being ohmic conductivity and the others related to the dielectric properties. 

We shall elaborate more about the modelling once we have examined various other 

problems related to the method. 

o* 

The discussion about c e ^ was centered around a specimen composed of two 
homogeneous materials. It can be easily shown that the results obtained are 
not altered if these materials are subdivided into slabs and stacked in series, 
while maintaining the condition that the total lengths dj and d 2 of the 
components remain unchanged. From the arguments leading to equation 05) it 
follows that if the first and last elements in the stacked sample are of the 
-same catena I , say I, the net charge stored in the sample is zero, whereas if 
they correspond to different materials the net charge, added up from all interfaces 
within the material, will be equal to the charge stored at the single interface. 



As a consequence of these considerations the arrangements in Figures 2(a) and 
(b) are equivalent as far as total impedance is concerned; the only difference 
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6 etween them is that in Figure 2(a) the net charge stored is zero whereas in 
Figure 2 ( b i s different from zero. 

Consider material 4 as made out of an arrangement of two homogeneous 

materials, I and 2, in contact (Figure 2c), while assuming that material 3 is 

* 

a good dielectric. The of material 4 will be given by equations (?*) D r 
i £3 ) , and accord i ng to the limiting case given by equation ( 2 fe) its Imaginary 

component will diverge when u)~*o . The question arises regarding possible 

* r* 

divergencies in the response of the composite material formed by £3 and £4 
since the last one presents the divergent behavior mentioned above. 

The total effective permittivity of the composite material in Figure 2(c) 
can be written as ( d 3 t d 4 ) £ 3 ^* 

r* 

'4 

* , + * * 

where £ 3 = £3 and C 4 . is the effective permittivity of and c z 

as given by equation (24). 

The limiting cases and ( 0 -*co show that no divergencies appear in the 

components of and £l{- . If materia! 4 was constituted of three or more com- 

C* 

ponents, an iteration of the process used to obtain c T shows that these results 
still hold. As long as the two end materials are good dielectrics the divergencies 
will be avoi ded . On the contrary, if material 3 is a metal (i.e. an electrode) 
a similar analysis (Alvarez, 1972, p. 29) shows that the divergent behavior re- 
appears when to-* o . 


£ r (% d 4 4 - £*<4 


The two domains of absorption 

In what follows material I will be the good dielectric and material 2 wi I I 
- represent the sample to be analyzed. In the usual treatments of the Maxwell- 

■4 

Wagner effect it is implicitly or explicitly assumed that the dielectrics i n- 

c" c * 

volved are non-di spersive (i.e. C t = C z = 0 ); this makes the conductivities 
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£>" ( and responsible for the Debye-type behavior, except for the 

divergence at low frequencies. If material I is a good dielectric the Debye- 

type behavior of the composite material is complete even at low frequencies 

(i.e. c e ^jr* o when co-*o ). This case has been treated by Arnoult et a I 

CI965). It is easily realized that the treatment of the Maxwel l-Wagner effect 

c" c" 

under the restriction c t 2 = 0 strongly limits its applicability. If this 
restriction is relaxed, only for one material, one can then analyze dispersive 
dreieelrics having non-negi i g i Die onmic conductivities such as rock samples 
and a series of heterogeneous semiconductors: ceramics, ferrites and the like. 

The only attempt, known to this author, to treat the Maxwel I -Wagner effect 
considering a dispersive dielectric is that of Dansas et al (1967). They pro- 
posed that if two materials of thicknesses d, and d 2 are placed together forming 

c* c / 

a sample, one of them being a good dielectric (i.e. C, = ) and the other 

presenting a Debve relaxation and ohmic conductivity ; i.e. 

£d2. - • Q 
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then, the effective dielectric permittivity of the composite sample can be 
expressed as 


£* - £■' + 
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-where sobTndexes d. and j3 refer to two different domains of absorption (Figure 3). 

Theo<-domarn corresponds to a low-frequency energy loss due to the motion 
of free charge carriers (i.e. due to the Maxwel I -Wagner effect) and the ^3 -domain 
corresponds to dielectric losses arising from the dispersion of material 2. 

They claim that the two domains plot as semicircles as shown in Figure 3. 
However, we have shown that this is an erroneous result (Alvarez, 1972, p. 46) 
-sfnce^ the point at which theoc-and ^-domains meet is not a zero; at most it is 
a minimum. Consequently the two domains of absorption cannot be represented by 


semicircles. 
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Instead of their expression for 8^ we obtained its components as 
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(2,6) 




by- subst ituti ng the given values of £* and £* in equation (24) ^ and with Q = ^ - 

Figures 4, 5 and 6 show plots in the comD lex p lane of these functions 
when specific values are assigned to the parameters involved. The existence 

of the two absorption domains is readily discernible. The data points do 

not plot in semicircles (dashed lines). 


Substances with a distribution of relaxation times 

In the previous discussion we have dealt with a substance presenting a 

Debye relaxation and a given ohmic conductivity; this is the case treated by 

Dansas et al (1967). This, of course, is a particular model and one should 

not be surprised of finding substances that do not behave in the way prescribed 

c * 

by the corresponding (i.e. equations (2ft) and U^)). But nothing in our 

previous treatment prevents us from generalizing the considerations made for 
the Debye substance to other types of relaxations. 

This generalization is of importance since it will permit us to treat sub- 
stances that present a distribution of relaxation times, rather than a single 
one as in a Debye substance. Mul ti-crysta 1 1 i ne materials, such as rocks, are 
likely to present this behavior. As illustrations of applicability of the 
method we have chosen the functions known as the Cole-Cole arc (Cole and Cole, 
1941 5 and the Cole-Davidson arc (Davidson and Cole, 1951) which represent cases 
of substances exhibiting a continuous distribution of Debye relaxation times. 
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The Cole-Cole arc is represented by 

c* = t ' + L ° ~ — o ^ i 

and Its plot corresponds to a semicircle whose center lies below the K^axis, 

the parameter h is a measure of how much the behavior differs from the 

correspond i ng Debye relaxation process (i.e. when h = 0). The Cole-Cole arc 

corresponds to a symmetrical distribution of relaxation times. 

Assuming that materia! 2 behoves in this manner and has ohmic conductivity 

G ^ , and that material I is a good dielectric we want to find the corresponding 
£ (i.e. equation 24). 
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are the components of when the dispersive sample presents a Cole-Cole 

distribution of relaxation times. 


The Cole-Davi dson arc is represented by 

* ' */ 

£ = * 4 ' (t+ju,z.yy 


o( 4z. I 


when plotted in the complex plane this function is a skewed arc, O. being a 

measure of skewness. This distribution of relaxation times is not symmetrical, 

rather it consists of contributions that diminish in importance as the frequency 

c" 

increases producing values of c smaller than the corresponding Debye relaxation 

(i.e. o( = I) towards the high-frequency end. Just as in the previous case we 

c* c* 

can obtain q ^ for the comDOsite substance in Fi-gure 7 when is described 


by a Cole-Davi dson arc; in this case 
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are the components of £ when the dispersive sample presents a Co I e-Davidson 
distribution of relaxation times. 

In equations (BO), (31), (32), and (33) we have shown that the generalization 
of this method of measurement and analysis may be carried out straightforwardly 
to substances presenting a distribution of relaxation times. 


The parameter 0 
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A discussion on the role of the parameter Q, the ratio of thicknesses 
^ x /^i, is necessary since it controls the magnitude of and . 

When analyzing the limiting cases of the general equation for 
we obtained equation (26) as 

Qdj +■ <3*i d. 

If material I is a good dielectric Q = 0 and the above equation- can be 




written as 


C* = () + §)£', £& — O . 

‘ff 


Thus, the maximum value of depends on Q and the (constant) value of the 

dielectric permittivity of material I. Once the latter is chosen, and 

a/ 

therefore its c, is fixed, it will be sufficient to vary its thickness through- 
out a series of experiments to obtain larger or smaller values of the static 
, while maintaining fixed the sample thickness. Obviously not only 
at zero frequency is affected by such a change in Q, but all other values 
°f SgQ and will change proportionally. The net result of having large Q-values 
(l.e. Q rj 20 or higher) will be the enhancement of the CA -domain while small 

values of Q (i.e. Q < 10) will permit a clearer definition of the ^-domain. 

c R 

Large Q-values imply larger values of both and C e ^ than the corresponding 

ones for low 0- values. 

_/ pH 

Thus the overall effect of Q on the values of c and CW is that: (I) 

when Q has large values one obtains magnified values of the components of C 
while an enhancement of the # -domain occurs, and (2) when Q has small values 
the components of present diminished values while enhancement of the 

ft -domain occurs. The meanings attached to magnified and diminished are, of 


course, relative; a more precise signification evolves when dealing with a 
Dartfcular set of parameters. 
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Varying the parameter Q implies varying the magnitude of the applied 


electric field in material 2 (i.e. the sample). This is the physical phenomenon 
responsible for the modification of the c e y values when Q varies (Alvarez, 

1972, pp. 52-56). 

For example, when decreasing Q (increasing d. and decreasing d y ) a smaller 

... *' hake 

electric field is established in material 2, cunfoi wi I ^decrease the losses in 
it (ohmic and dielectric) as well as in the composite sample. The response of 
material I is not affected since it is a lossless dielectric. 

The experimental ability to modify the response of a given sample without 
altering its geometry or electrical properties, as well as the possibility of 
enhancing a given domain of absorption represent a definite advantage to the 
study of rock samples since they allow us to make a more detailed study of the 
regions of interest, in particular the low-frequency region which is the one 

that traditionally has presented more experimental problems. 

4k 

Effects of parameters variation on Key plots 

A few examples shall illustrate the variations of due to variations 

of the parameters involved. Relative dielectric permittivities ( K*. e /e - > 

will be used since they are more commonly found in reporting experimental results. 

Six or seven parameters will characterize each example: IC] will be the 

(constant) value of the relative permittivity of material I; K' z and K w i will 

be the extreme values of the K of material 2; Z> 2 the dielectric relaxation time 

of material 2 and its ohmic conductivity; Q is the thickness ratio d^/dj. A 

Debye substance is completely characterized through equations (28) and (2^) and 

the parameters listed above. If a Cole-Cole arc or a Cole-Davi dson arc are 

•modeled, an additional parameter is necessary: h for the former and CX for the 

latter. The corresponding equations are (30} , (31), (32) and (33). 

* 

Table I establishes the relation between the figure number and the 
corresponding parameter values. In Figures 4 through 13 the frequencies are given 

in kilohertz and the axes are simply labeled % and K. instead of anc * ^ e ff 
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Figures 4, 5 and 6 constitute an example of how the parameter Q affects 
'K 

the piots of in the complex plane when materia! 2 is assumed to behave as a 

Debye substance. Q takes up values 4, 12, and 48, white the rest of the para- 
meters are J<ept fixed at values: K t = 6.0, z = 4.0, = 250, = 7.95 x I0~ ^ 

sec, and (^ 2 = 1 .0 x 10 ^ *. The value K' 2 = 250 may seem too high if 

one thinks in terms of a homogeneous substance but when substance, 2 is con- 
sidered heterogeneous the Maxwe I l-Wagner effect will give rise to that kind of 
values for the static dieiectric permittivity. The value of Z 2 = 7.95 x !0 ^ 
sec corresponds to a critical frequency of 200 kilohertz. 

The presence of thecx- and domains is clearly discernible in Figures 
4^ 5 and 6. The frequency range 30 Hz- 100 MHz was chosen for the theoretical 
examoles because it corresponds to an easily accessible range with commercial 
bridges. Notice how the o( -domain becomes enhanced in going from low to high 

values of Q. The dashed lines correspond to the semicircles proposed by Dansas 

* 

et at (1967); obviously the K ^ plots are not semicircles. The point at which 

„ " IS* 

the £X- and -domains meet presents a non-zero, although small, K ^ value. 

Figures 7 and 8 show the effect of variation in the conductivity ^ of a 

sample presenting a Cole-Cple distribution of relaxation times (equations 30 and 

3M. A tendency is observed to fuse the two domains of absorption into only one 

when the ohmic conductivity increases; this trend is observed in other types of 

-substances (i.e. Debye, Cole-Davi dson, etc.) whenever Q increases. A shifting 

of the points towards higher values also takes place. 

» 

Figures 9 and 10 show the effect of variation of the parameter o( in a Cole- 
Davidson distribution (equations 32 and 33)- The closer is <X to unity the closer 
will be the behavior of the Col e-Davidson distribution to a Debye relaxation. 
Theo<-domain remains practically unchanged upon o< variations while the ^-domain 
is drastically modified. 



A comparison can be made between Debye and Cole-Davi dson relaxations 
with the aid of Figures 5, 9 and 10; all parameters, exceptC*, are the same 
in these figures. A more complete evaluation of the effects of parameters 
variation has been made elsewhere (Alvarez, 1972, pp. 56-64). 

EXPERIMENTAL DISCUSSION 

We want to show now an experimental instance in which the two domains of 
absorption appear clearly d< fferen+i a+ed- Three sets of measurements will “'be 
reported, all of them made on a sample of hematitic sandstone at room conditions. 
Each set corresponds to a different Q-value so that the predicted variation in 
the response of the sample can be checked experimentally. 

The only aim we have in mind in presenting these examples is to exhibit 

the applicability to rock samples of the method discussed herein. Future work 

designed specifically to obtain handbook- type values of various rock specimens 

will have to include a complete mineralogical description of the samples as 
wefl as full information on the particular experimental conditions under which 

the measurements were made. Such studies lie beyond the scope of the present 

work and will not be found In the examples to be presented. 

Two Impedance bridges were used in the experimental measurements; a General 
Radio CapacTtance Measuring Assembly Type I6I0-B was used to analyze the 30 
Hz-fOQ KHz range and a Boonton Q-Meter Type 160-A was used to analyze the 100 KHz 
20 MHz range of frequencies. The former has an accuracy of +0.1 % _+ 0.6^ for 
capacitance measurements and a dissipation factor (D) accuracy of 0.00005, 
while the latter has +_ 2% and +_ 10^ nominal accuracy for capacitance and Q-measure 
ments respectively. Further descriptions of these instruments can be found in 
the manufacturer’s literature (General Radio, 1965, and Boonton Radio). 

•The General Radio Assembly includes a type 7I6-P4 guard circuit which 
permits three-terminal measurements to be made. All our measurements in the 
range of frequencies covered by this bridge were made in the three-termina I sub- 
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stltution method of measurement; therefore all contributions arising from 
parasitic capacitances were cancelled out and the measurements in the 30 Hz - 
100 KHz range are considered free, within the error limits stated, of such 
contributions. The sample was located, in all measurements. In an electrically 
shielded enclosure to avoid noise pick-up. 

The Boonton Q-Meter does not have capabilities for three-terminal measure- 
ments and consequently there is the possibility of contributions from parasitic 
capacitances, especial ly at the higher frequencies. 

Considering all sources of error. Instrumental and geometrical, we have 
estimated the overall accuracy in our measurements as _+ 5% in the 30 Hz - 100 KHz 


range for and Ky ; 10$ for /^and +_ 15$ for K^in the 100 KHz - 10 MHz range, 

and +_ 15$ and + 20$ for /^and /^.respectively in the 10 - 20 MHz range. 

As the good dielectric we used mica sheets of various thicknesses. The 
metallic electrodes were made of Indi um-Mercury (In-Hg) amalgam which has proved 
to be an excellent substance for making intimate contacts with rock samples in 
general although in the present cases it was placed in one surface of the mica 
sheets . 

The thicknesses of the mica sheets and sample were independently determined 
with a micrometer of 1/100 mm/div accuracy. Several measurements were made at 
various places for each dimension and then these were averaged to obtain a 
representative value. 

The values for the dielectric constant of mica ( Kj ) were experimentally 
determined, over the range of frequencies studied, prior to mounting them in the 
composite sample. At the low frequency end (i.e f 300 Hz) the losses in the 
mica usually increased by about an order of magnitude; typically A. =0.18 at 
30 Hz which in all cases was the lowest frequency employed. Three different values 
were'obtai ned for the three cases analyzed: = 4.11, 6.70 and ^ = 6.90. 

It should be mentioned that the mica samples used came from at least two different 
types of mica; one of them presenting black streaks while the other was clearer 


and of a uniform color. 
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Table 2 displays the geometrical parameters involved; the rock sample 
was the same in all three cases; d| is the mica thickness, 62 the sample thick- 
ness, Q is the ratio anc * A is the area of the electrodes. They are re- 

lated to the figure number In which the experimental results are found as 
crosses joined by a thick line; the numbers by the crosses represent the fre- 
quencies, in kilohertz, at which the value was obtained (Figures II, 12 and 13). 



experimental curve is maintained. The point at which ihec\-and ^-domains meet 
is around f = 0.8 KHz. 80 th domains are only partially sketched by the 
experimental points, due to experimental limitations at both ends of the fre- 
quency range, but there is no doubt about their presence. 

A very interesting feature appears in Figures 12 and 13; this is a strong 
and narrow absorption which was missed in the first experimental measurements 
with Q = 4.82 when we went from 700 KHz to 2 MHz; notice that it is only inferred 
in Figure II, but in the next two experiments (i.e. Q = 9.54 and Q = 24.1) the 
peak appeared quite clearly. In fact we had to repeat the experiments correspond- 
ing to these two last values of Q, measuring again around I MHz, in order to make 
sure that it was not a resonance due to some particular geometrical configuration 
of sample holder and leads, as was pointed out to us as a possibility by 
Westphal (1971). ^ '"Wfe hold the view that the peak is a 

true absorption peak associated with the hematitic sandstone sample studied, 
not a false absorption due to instrumentation, for the following reasons: (?) 

the peak appears only for hematitic sandstone and not for other samples we have 
stydied f although not reported in this work, 

and (ii) the peak is reproducible not only when repositioning the sample in the 
Sample holder but when the overall thickness of the composite sample is varied. 
Since any explanation of the presence of this absorption peak would have to be 
made in terms of a particular group of molecules present in the sample we shall 
-offer no such explanation at this point. 



In Figures II, 12 and 13 we have included, in addition to the 79 

experimental data, the theoretical points corresponding the the same 
frequencies used in the experiments as obtained from equations (28) and (29). 

We assumed that the experimental response was due to a sample in which there 
was a distribution of relaxation times dominated, in the low frequency region 
(i.e. f < 500 KHz), by a Debye process. This was suggested by the shape of the 
experimental curves. Since the high absorption peak occurs at frequencies around 
hMHz,. and since there are only a few experimental points above it, we con- 
sidered it unnecessary for the purpose of these examples to try to fit another 
theoretical function to the experimental data. 

The theoretical points are shown as empty circles joined by a thin line 
and- the numbers correspond to the frequencies. The theoretical parameters AT, 
and' Q are the same, as in the experimental cases; they were kept fixed while 
5, , and were varied, using the functions given by equations 
(?&) and (2.9), to fit the experimental points. The fit has not been optimized 
SFincerour purpose is only to show the mechanism by which the parameters corres- 
prandi ng r to the rock sample can be extracted from the experimental data. Once 
the" set" of parameters = 80, ^oz~ 400, = 4.0 x 10 ^ sec, and = 4.52 x 

1 0T. (Q ;W) * was fixed for Q = 9.54 and the theoretical points in Figure 12 
attained, the values of Q and Ay were changed to those shown in Figures II and 13, 
si that" the corresDondi ng theoretical curves could be obtained. Thus in Figures 
II, IZ and 13, i-f is shown how well a single set of theoretical parameters 


correspond simultaneously to the three sets of experimental data. Notice that 
the theoretical minimum at K - correspond closely to the experimental 

minima, although the former tends to concentrate several frequencies around such 


paint" while the latter shows a smoother spreading of them. 

The. parameters K 0l and ^ 2 had to be assigned large values in order for the 
theafe+ical points to fall at the positions shown; smaller values (i.e. 

K Ai = 200 and j ^ = 40) gave theoretical point positions not related at all to the 
experimental ones. The value of 400 for ^, 2 is attributed to the Maxwel l-Wagner 
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effect taking place within the rock sample and should be contrasted to the 

3 

values of K.a/|o at I KHz obtained when we made direct contacts between 
electrodes and sample. The value K z = 80 suggests the presence of a water 
effect in the sample, but this constitutes only a speculation at the present 
stage of analysis. 

Although the re ( are discrepancies that indicafe that the assumed Debye 
behavior at the low frequency end is at best an approximation to the real 
oenavior ot the sample in Such a region (as the one just discussed for the 
distribution of points around ), such discrepancies tend to be the 

same in the three cases studied. 

We want to emphasize once more that the conductivity 0^ obtained with this 
method refers only to the ohmic conductivity of the sample. The contribution 
to the losses arising from the dielectric properties, that is, the dielectric 
conductivity, can be obtained for every frequency from the parameters , 

KLx and and the adequate functions describing its behavior U.e. Debye, 
Cole-Cole, etc.). 

We have presented, in this experimental discussion, evidence of the 
applicability of the method proposed in this work to the measurement of rock 
samples. Also we have tried to suggest some possible ways of analyzing the 
response of the rock, but obviously much more research is needed in this respect 
to fully exploit the possibilities furnished by the method. 


CONCLUSIONS 

The similarity in electrical behavior between rock samples and heterogeneous 
dielectrics with non-negl igible ohmic conductivities has led us to a general 
analysis of the Maxwel l-Wagner effect; from this it has been concluded that in 
— order to have electrical steady state conditions in that type of materials it is 
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necessary to have volume charge accumulations at the interfaces within the 

... 

sample. The expression for the charge density ( 2. ) accumulated at an inter- 

face as a function of frequency has been obtained and shown to plot as a semicircle 

In the complex plane. Expressions for the effective dielectric permittivity 

c* 

( Cfff) a composite material have been presented showing that the relaxation 

— * 

time associated to it was the same as the one for the charge density 2. when 
the sample components were considered to be non-dispersive dielectrics of differenl 
Ohmic conductivities. 

¥ 

The analysis of the divergence of the imaginary component of 8^ has made 
it clear that such divergence arises whenever the two components of a composite 

Sample have non-negl igible ohmic conductivities, whereas if this parameter 

is negligible for one of them the divergence at low frequencies is;avoided. 

This has made it possible to develop the method of measurement for heterogeneous 

dielectrics having non-neg I i gi b le ohmic conductivities by requiring that one 

of the two components in a composite material be a good dielectric. 

The Ideas advanced by Dansas et al (1967) about the behavior of the function 

were discussed. We indicated that their conclusion about the two domains 

of absorption plotting as semicircles in the complex plane does not hold; 

instead the correct expressions have been obtained. Furthermore the method has 

been generalized to include dielectric substances presenting a distribution of 

relaxation times instead of the single one involved in the Debye behavior. A 

few examples on three materials presenting different behavior (Debye, Cole-Cole 

and Co I e-Davidson) have been analyzed to illustrate the type of responses to be 

expected from this method. Three experimental examples have shown the feasibility 

of studying electrical rock properties by means of this method. 

As can be appreciated, the method presents the advantage of permitting’ 

application of the pertinent theory already existent for. dielectrics; we have 

only touched upon this matter when extending the di scussion to substances having 



a distribution of relaxation times. There are many other interesting possibilities 
such as that of evaluating the effects of a distribution of ohmic conductivities 
In the material . 

The need for further work in the area of interpretation of results as 
obtained with this method has been stressed; meanwhile the possibility of 
obtaining new data on rock samples is opened, as may be the case with many 
sulfides (i.e. Fe$2> PbS, C^S, etc.) whose dielectric properties are almost 
unknown and poorly described. The main reason for such a lack of data is 
their relatively high conductivities which heretofore have prevented their 
measurement in the frequency region of interest in exploration geophysics. 
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TABLE I. Parameters used in the computation of examples in which 

the dispersive dielectric is characterized by Debye, Cole- 
Cole and Cole-Davidson relaxations. 
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rig. 
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(mm) 

^2 

(mm) 

A 2 
(crn ) 

Q 

II 

0.85 

4.10 

3.84 

4.82 

12 

0.43 

4.10 

3.84 

9.54 

13 

0.17 

4. 10 

3.84 

24. 10 


TABLE 2: Dimensions of sample electrodes and 

• micas used in the experimental examples. 
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Fig. I 


Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Fig. 6. 


Fig. 7. 


Response of the total current density (i.e. conduction plus displace- 
ment currents) in a sample formed by two slabs of different materials 
when the voltage shown is applied. 

(a) The net charge stored in the two interfaces in the steady state is 
zero, (b) The sample has materials 3 and 4 lumped together, there 
Is only one interface and the net charge stored in it is different 
from zero. The total impedance remains unchanged, (c) Material 4 

is formed by two homogeneous components characterized by £. 

c< c * 

and c x ; represents a good (i.e. non-lossy) dielectric. 

Cj^plot in the complex p I ane, accord i ng to Kansas et al (1967), of 
of a sample formed by a lossless dielectric and a material presenting 
a Debye relaxation and non-zero ohmic conductivity, showing the two 
domains of absorption as semicircles. The relative sizes of the 
°C - and 0 -domains can be varied by changing the thickness ratio 

d*/d, 

K, variations arising from changes in Q. Compare to Figures 5 
and 6. Material 2 is described by a Debye function. The dashed 
fines correspond to semicircles. 

K -variations arising from changes in Q. Compare to Figures 4 

TJ 

, I 

and 6. Material 2 is described by a Debye function. The dashed 
lines correspond to semicircles. 

^variations arising from changes in Q. Compare to Figures 4 

D 

and 5. Material 2 is described by a Debye function. The dashed 
lines correspond to semicircles. 

/(^variations arising from changes in . Compare to Figure 8. 
'Material 2 is described by a^Cole-Cole function. 

/^variations arising from changes In • Compare to Figure 7. 
Material 2 is described by a Cole-Cole function. For the same 


Fig. 8. 
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frequency, a point acquires larger values when the conductivity 
increases. 

Fig. 9. ^variations arising from changes in the parameter ck . Compare 

to Figure 10. Material 2 is described by a Coie-Davi dson relaxation. 
Compare to Figure 5 described by a Debye relaxation. 

Fig. 10. K variations arising from changes in the parameter (X . Compare 

n 

to Figure 9. Theo<-domain remains unaffected while the ^-domain under- 

o 

does a cons • de r ab *e f ir»+inn. Material 2 is described bv a Cole- 
Davidson relaxation. Compare to Figure 5 described by a Debye 
re I axat i on . 

4 

Fig. II. Experimental values of for a hematitic sandstone sample, of 

2 

3.84 cm area and 4.10 mm thickness, are indicated by crosses joined 

by the thick line. The points shown as empty circles correspond to 
a Debye substance characterized by the theoretical parameters in- 
dicated; these are the same for the three experimental cases and 
therefore show how well a single set of theoretical parameters corres- 
pond simultaneously to the three sets of experimental data (see 
Figures 12 and 13). It was assumed that the experimenlal plot was 


Fig. 12. 


dominated by a Debye relaxation in the frequency region £ £00 KHz. 
Experimental values of ^ for the same sample as in the previous 


in n- et 


figure. The value of Q is 9.54. There is an increase in ^ ef/ 

1 / ^ 

and values with respect to the plot of Q = 4.82. A strong 

absorption peak is delineated by five experimenlal values around 1.4 MHz; 
in the case Q = 4.82 this absorption was not detected due to the lack 
of measurements between 700 KHz and 2MHz. Its presence has only been 
inferred in that case. 


* 
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♦ 

Fig. 13. Experimental values of Kg,, for the same sample as in Figures II and 



from the previous figures. The strong absorption peak is delineated 
by four points around 900 KHz. The experimental curve is smoother 
than the previous ones; this effect is thought to be due to a 
different atmospheric moisture content at the time of measurement. 
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Electrical Non-linear Behavior In Rocks 
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ABSTRACT 

Previous studies on rocks saturated with electrolytes have shown that 
there is an upper limit for current density, up to which linear behavior is 
observed. 

Evidence of the opposite behavior in some dry resistive rocks in contact 
with metallic electrodes is presented, showing the existence of a non-linear 
region followed by a linear one when the voltage, or current density, is 
Increased. 

This effect is attributed to interfacial phenomena occurring at electrode- 
sample interfaces. It has been shown elsewhere that metal-insulator or metal- 
semiconductor interfaces may give rise to linear or non-linear behavior 
depending on the type of electric contact established: ohmic, blocking, or neutral. 



Five resistive and two conductive samples were analyzed. Of the resistive 
group, limestone, granite and hematitic sandstone showed a non-linear region, 

whereas argillaceous sandstone and muscovite schist did not, suggesting that 

; 

different type of contacts are extablished by the same metallic electrode 
Cindium-mercury amalgam) in the various samples. Galena and pyrite, 
samples of the conductive group, behaved iinearly f as expected. 

The existence of small polarization voitages was observed in the resistive 
group of samples; ttwy uiere attributed to space charge concentrations at the 
electrode sample boundary. Experimental determination of the potential 
distribution on a sample surface reveals high non-uniformity of the electric 
field in that region as well as the presence of space charges immediately 
adjacent to the electrodes. 

The polarization voltages observed give rise to a displacement of the 
current vs voltage characteristics from the origin. It is shown that if 
the I -V characteristics do not go through the origin, the corresponding 
resistivity vs current density plot will appear to be non-linear even 
when the l-V characteristic is linear. A warning is made regarding the use 
of resistivity vs current density plots to decide on non-linear behavior. 

Finally, it is suggested that interfacial phenomena and non-linear 
effects may be common features in two- and four-electrode systems. 



Electrical Non-linear Behavior in Rocks 


INTRODUCTION 

In the study of electrical properties of rocks the assumption of 
linear Behavior is often made, in addition, it is admitted that rocks exhibit 
linearity up to a given current density value; above it non-linear effects 
.<noear. The existence of an upper limit in current density has been 
established in studies such as those of Madden and Marshall (1959) and 
Anderson and Keller (1964). 

information on the subject of non-linearity is scant and often 
incidental to other studies of electrical conduction in rocks. McEuen et a I 
(1959) observed, for synthetic raeta I I if erous ore, an impedance that was 
deoendent on current density. Scott and West Cl 969) observed a decrease of 
resistivity with increasing current density in samples containing disseminated 
ssffides. Zonge (personal communication) has observed non-linear effects 
W-.ile studying IP responses, and Katsube and Collett (1972) have recently 
reoorted on electrical non-linearity using ac excitation. Possibly the only 
s/stematic approach to the problem has been that of Shaub (1965 and 1969; 

^•'aub and Ivanov, 1971). 

While studies in electrochemistry (.Shaw: and Reraick, 1950) may 
^•/e as a guideline to understand non-linear effects in rocks in which 
y -eous electrolytic conduction predominates, care should be exercised not 
"''Extrapolate results to samples in which such conditions may not be applicabl 

When sufficiently intense electric fieldsare applied to a material, 

1 -rill eventually undergo an irreversible change in electrical properties; 

>re arriving at such a limit electrical non-linear phenomena will appear. 



108 


A system described electrically by the relation 

J- o-L, 

Cl) 

in which 6Wtf(E), exhibits a linear (or ohmic) behavior. J is the 

9 

current density, the electric field and 0^ the conductivity.* 

If the conductivity O’ , in equation (I), is a function of the 
electric field, tnat is, if g"— (tCE), the relation between the current 
density and the electric field can be generally written as a polynomial 

yj = (f'E E •f-d'"' 11 E +- 

II (2) 

in which the coefficients <K (F , ..are constants. Relation (2) 

describes a non-linear behavior. 

If the total electric field is the sum of two independent 
electric fields and E^i r + Is readily verified that superposition holds for 
relation ( I ) / _ r \ 

J- <Kf. *Mx) 

whereas for equation (2) superposition no longer holds: 

j . ( e, + f,) + J 1 " ( f , 2 * d * &A) <■•••• 


Consequently, if the fields E_| and Eg are of the form 

ICO, t J U)q,E 

k>! ~ and £z ~ £ 


(3) 


harmonic generation, in equation C3), will occur (Shaub, 1965). 

If relation (2) describes the behavior of a given material, then 
under small electric fields, and if the condition 

G~'yy 6-"» <p! H 

* MKS units are used in this paper. 
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holds, the response <wi 1 1 be linear. Clearly, for sufficiently large E_ fiel ds, 

0-"E 2 ^ fl-'E 

and the non-l inearities will appear. 

However, intense electric fields are not the only possible cause 
of non-l inearities. Non-linear behavior may also be observed at small elec- 
tric fields in insulators and semiconductors in contact with metallic 
electrodes. 

Contact effects have been discussed in relation to insulators and 
semi conductors by Mott and Gurney (1950), Frank and. Simmons (1967) and others. 
Simmons (1971) has reviewed the three types of contacts that may arise in 
metal -insulator-metal and meta l-semiconductor-meta I systems: ohmic, blocking 
and neutral. Once an ohmic or blocking contact is established, linear and 
non-linear behavior may be obtained for the same meta I -insulator-metal system 
depending on the voltage applied. Furthermore, a contact that in a given 
voltage region is ohmic may, at different voltages, become blocking and 
viceversa. A space charge is often found at the meta I -insulator or meta I - 
semiconductor interface; it plays a central role in the contact properties. 

Consequently, the possibility exists of observing similar phenomena 
in the case of rock specimens, which are either semiconductors or insulators. 

In contact with metallic electrodes. 

Anderson and Keller (1964) established the existence of an upper 
limit for current density, up to which the behavior of a pyrite sample was 
linear. They used direct current in this determination. We have also used 
direct current in the experiments reported herein since we are looking for 
electrode effects as the origin of non-linearities ; they will be most clearly 
shown under these conditions. 

The measurements were performed with samples exposed to room 



conditions although, for comparison purposes, a few were made at high, vacuum 

conditions (i.e. I0~^ Torr). At room conditions the samples were , of course, 

subject to the effects of atmospheric mo isture, whTch ranged from 20 to 45iS 

relative humidity. A two-electrode system was used throughout the experiments. 

arcftlliiceous zanJs bane., 

The samples studied can be divided in two groups:^ granite, 

1 1 roes tone and heraatitic sandstone belong to the group hereafter called 
"resistive samples", and samples whose main constituents are galena and pyrite 
respectively belong to the "conductive samples". 

NON-LINEARITIES: EXPERIMENTAL 


The dc regulated power supply used was a Princeton Applied Research 

model TC-602CR. Total current and voltage in the sample were measured with 

14 

two electrometers (Keithley, models 601 and 6I0B) of 10 ohms input impedance. 
A two electrode system was used; Indium-Mercury electrodes were directly 
applied to the sample in each case and, as a result, a Maxwel l-Wagner effect 
was observed when measuring the samples belonging to the resistive group 
(.Alvarez, 1972). Transients occur when establishing or discontinuing dc 


excitation; the measurements reported were made after all transients had 


•i 


disappeared. 

It has become common practice to report values of resistivity ^ , 
for rocks, derived from measurements of total current and voltage in the sample 
When the sample geometry corresponds to parallel faces of area A and thickness 


d the relations 


V 


Y. - R = P ^T* 

l ~ i A 


are used to obtain ^ . This operation results in deriving a point prop- 

erty ti.e. ^(x, Y* z) from the bulk properties of the sample (i.e. total I 
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and V). Clearly, this way of proceeding will be valid only if the field E 
in the sample is uniform. Further ahead w;e shall present experimental evidence 
showing the non-uniformity of the field on the surface of a sample. Consequently, 

i 

we have tried to avoid the above operation wherever possible. The data 
is presented mainly as l-V or R-V plots ( Figures la and b); the measured 
parameters were current and voltage in all cases. Since the sample geometries 
are given, one may obtain, if necessary, the sample resistivities from the 
resistance plots by making the assumption of uniform field in the sample; at 
any rate, they will differ only by the geometric factor d/A. Two plots of 

resistivity vs current density are included to a I low compari son with previous 

A 

work (i.e. Anderson and Keller, 1964, and Scott and West, 1969). 

In Figure 2a the total resistance R of a limestone sample of the 


dimensions shown is plotted against the total 


^voltage across the 


sample. The shaded surface represents an electrode. For voltages of less 
than 5 volts a non-linear region is observed; linear behavior appears to 
occur for larger voltages. Notice that for V = 0. I volts a negative value 
of R is obtained. This was due to current flowing in the direction opposite 
to that determined by the polarity of the applied voltage. The only way such 
an effect can be accounted for is by the existence of a polarization within 
the electrode-rock-electrode system , larger than and of opposite polarity 
to the one externally applied. 

In Figure 2b data corresponding to a similar sample of limestone 
is plotted. The voltage goes from -8.0 to + 7.0 volts; in this case one 
electrode of the sample was constantly grounded while the other was varied 


from positive to negative voltages. The presence of a polarization within 
the electrode-rock-electrode system is again manifested by a zero current 


flow a.t a voltage of +0.18 volts. 


Correspondingly, there is 



a non-zero current of -0.3 x 10 amp when the sample terminals are short - 
cTrcuited. 

The data in Figure 2&jis plotted as resistivity vs current density 
modulus in Figure 2c; according to the established criterion ( Anderson and 
Keller, 1964), the resistivity varies non-linearly with current density. A 

discussion of these results and the way of plotting them is made below. 

The fact that a voltage exists between the sample Terminals when 

no excitation is applied constitutes an important difference with respect to 
the common meter ia Is obeying Ohm's law; in these the current vs voltage char- 
acteristic is a straight line passing through the origin, whereas in the 
rock sample it is displaced from the origin (Figure 2b). 

In Figure 3a the current-voltage characteristic (the scales of current 
and voltage are multiplied by arbitrary factors) of a material assumed to have 
such properties is used to model the behavior of the limestone sample. They 
show a voltage of I unit when the terminals of the model sample are open- 
circuited, and a current of -2 units when they are short-circuited. Other 
than the peculiaritv of being displaced from the origin the characteristic 
is linear. Assuming that the geometry of such a sample is a unit cube 
we can obtain the resistivity vs current density plot (Figure 3b) correspondi ng 
to the characteristic in Figure 3a. 

In Figure 3b there are two branches of resistivity values, denoted 
P and § . Since resistivity is obtained by making the ratio V/l and 

multiplying by a geometric factor (in the present case I), it will be positive 
for voltages V>1 and V<0, in which both V and I have the same sign; for 
the region 0<V<l the voltage is positive, but the current is negative, thus 
the resistivity will be negative. It is for this reason that we have been 
forced into defining negative resistances and resistivities. In common mater- 
ials obeying Ohms' law this problem never arises. 



For the same absolute value of current, in Figure 3a, there are 


two different values of voltage. Therefore there will be two different 

values of resistivity for each value of the current density modulus 

CFigure 3b). This is true up to a value of current density, approximately 

100 times larger than the one corresponding to zero resistivity C i.e. to 

the current density J 0 flowing when the sample is short-circui ted. ) Above 

such a value of current density the resistivity corresponds to that of 

i-e. 

the reciprocal slope of the characteristic (^the actual resistivity). 

Thus, a material with a linear l-V characteristic, displaced 
from the origin, gives rise to an apparent non-linear behavior when plotted 
as resistivity vs current density. The apparent non-linearities appear 
for values of the current density modulus IJI CI00 J e ; whereas for |j| > 100 
the resistivity is linear. 

In trying to determine non-linear behavior with dc excitations, 
it is important to establish whether or not the l-V characteristic of the 
rock sample is displaced from the origin; otherwise non-linear behavior may 
be inferred from resistivity vs current density plots for samples with 
linear l-V responses. Therefore, it seems more adequate to decide on the 
linearity of the rock response on the basis of l-V rather than J-J plots. 

Returning to the results of the limestone sample one can observe 
(Figure 2 fc>) the displacement from the origin plus an actual non-linear plot 
of the l-V characteristic. Figure 2c shows the behavior predicted from 
Figures 3a and 3b; except that at the higher current densities the ^ and §> 
values, although coincident, show a non-zero slope; this fact plus a distor- 
tion of the curve shapes with respect to the theoretical one (Figure 3b) 
arise from the non-linear behavior of the limestone sample. 
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Figure 4a shows the l-V plot corresponding to granite. A linear 
region, for y >3 volts, and a non-linear region, for V4 3 volts, are 
obtained. Figure 4b, the corresponding R-V plot, suggests a non-linear region 
below 3 volts. Notice, however, that from Figure 4b alone one would not 
be able to discriminate between the increase in R, when V decreases, due to 
an actual nnn-ifnear effect or to that due to a displacement of the l-V 
characteristic from the origin. 

The main difference of the granite data with respect to the lime- 
stone data is that the latter shows non-linear behavior throughout the voltage 
region analyzed, whereas the former presents a linear region from 3 to 8 volts. 

The l-V plot corresponding to an argillaceous sandstone sample is 
presented in Figure 5. The characteristic is linear, and, in the scales 
plotted, appears not to be displaced from the origin. 

In Figure 6 data corresponding to a muscovite schist sample plots 
as a straight line; the behavior is linear but there is a small displacement 
of the characteristi c from the origin, towards the positive side of the 
current axis, contrasting with that of the limestone sample (Figure 2b), 
which intersected the current axis on t^ie negative side. As a result, the 
resistance values will decrease, instead of increase as in the limestone 
sample, when zero voltage is approached from the positive side of the voltage 
axis . 

Figures 7a and b correspond to the R-V plots of argillaceous 
sandstone and muscovite schist respectively. In Figure 7a a linear region 
seems to start at around one volt; for smaller voltages the decreasing 
R-values suggest a slight displacement from the origin of the l-V character- 
istic. This feature is masked in Figure 5 as a consequence of plotting 



in the scales shown, which are too large for the smaller values. 

The R-Y plot for the muscovite schist CFigure 7h) gives no clear 
indication of I inear Behavior, contrasting with the obvious linearity shown 
by the same data when plotted as an i-V characteristic (Figure 6). This 

W 

example underlies the importance of not relying on R-V or 0-J plots alone 
to analyze the linearity of rock samples. 

The points shown as crosses in a circle, in Figures 7a and b, 
were taken immediately after having applied the largest voltages to the 
samples (i.e. 9 volts for the limestone sample and 8 volts for the muscovite 
schist sample). In Figure 7a three such points delineate a behavior similar 
to the one observed in the case of limestone (Figure 2a), and in Figure 7b 
the corresponding point suggests a similar trend. This points out that a 
shifting of the l-V characteristics, from its original position, was brought 
about In each case by the application of the external voltage. 

The change is thought to consist of a redistribution of charge 
carriers within the material, caused by the passage of current. However, 
the change is not permanent since the sample presents the peculiarity of 
recovering its original properties if no new excitations are applied in 
several hours. Zheludev (1971) discusses the factors that may produce 
redistribution of charge carriers in a material and their relation to the 
space charge accumulation usually accompanying it. 

Figures 8a, b, c, and d correspond to a set of measurements on a 
hematitic sandstone sample. Figure 8a represents the typical response of a 
"resistive” rock sample, measured with a two-electrode system, when an external 
voltage V 0 that has been applied to the sample for some time is interrupted. 

The response consists of a "fast drop" in voltage (Vp) followed by a voltage 



decay resembling an exponential function. This remaining voltage (.y^) 
decays to a small value Cwith. respect to the magnitude of the applied 
excitation) within a few seconds, but seldom goes to zero, even after 
several hours of observation; it seems to correspond to the displacement 
from the origin of the i-V characteristics. 

The measurements on the hematitic sandstone sample include a 
comparison of the total resistance of the sample, designated Rp (obtained 
as the ratio of the total applied voltage V q to the total current Ip in 
the sample; Figure 8c), to the resistance Rp (Figure. 8b) defined as the 
ratio of the voltage Vp (Figure 8a) to the total current Ip in the sample. 

As can be appreciated from Figure 8b, there is a linear range for 
the resistance corresponding to the fast decay (Rp), starting between 3.0 
volts and 4.0 volts, whereas the total resistance (Rp) in Figure 8c, shows 
a non-linear behavior up to 9.0 volts. Below 3.0 volts Rp becomes non-linear; 
the negative value of Rp at V=0.5 volts requires additional description: when 
the excitation was interrupted, an increase in voltage took place 

Ci .e. instead of the usual fast decaying voltage, there was a fast increment) 
it was defined as negative Vp. This increase was followed by the usual 
decay. The measurement was repeated at least four times and the response was 
reproduci ble.^ Figure 8d shows the l-V characteristic obtained for the 
hematitic sandstone. Non-linearity is observed below 3 volts; above 3 volts 
the behavior is linear. The l-V characteristic shows a displacement of +1.0 
volt from the origin. On the basis of further verifications, a value of +0.5 
volts seems more adequate than the +1.0 volt value obtained in these 
measurements. The value of +0.5 volts was obtained after periods of several 
weeks of not subjecting the sample to voltage excitations. Observation 
times ranging from a few minutes to half an hour were necessary to insure 
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that steady state was attained when obtaining the data of Figures b, c, and d. 

A 

It appears that this hematitic sandstone sample requires much longer periods, 
than the other resistive samples analyzed, to return to its normal state 
after voltage excitations have been applied to it. Periods of around one 
hour necessary to attain equilibrium conditions have also been reported in 
ac measurements (Scott and West, 1969). 


A discussion of a schematic model shall be made in relation to 


the results of Figures 8a, b, c, and d, which may help to understand the 
origin of the non-linearities. One may think of the two decay regions in 
the voltage response as due to two different relaxation mechanisms, 
corresponding to at least two different species of charge carriers having 
a marked difference in relaxation times. The fast decay is similar to 
that of samples in which electronic conduction predominates, whereas the 
slow decay resembles ionic diffusion characterized by low mobilities. 

Although the latter mechanism may account for displacements of the 

2 

I — V characteristic arising from the application of external excitations , it 
cannot account for the di so lacements observed when the sample has not been 
subjected to voltage excitations at all, or when the sample is not excited in 
periods of several days: these need to be explained in terms of a mechanism of 
different origin. 

Information on the same type of hematitic sandstone reported elsewhere 
(Alvarez, 1971) shows that the magnitude of the fast voltage decay (V ) , attributed 
to electronic conduction in this model, tends to zero when the atmospheric moisture 
content is reduced, by evacuation, in the sample; eventually, in a vacuum of 
around 10 ^ Torr, V_ becomes zero. At this time the decay voltage V D equals the 
applied voltage V 0 , since-Vp = 0. 


2 


See explanation related to the crossed circles in Figures 7a and b. 
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In terms of the model under discussion, Vp = 0 implies that the electronic 
conduction mechanism has been supressed, leaving only the low-mobility charge 
carriers to conduct current through the sample. But since the water has been 
removed from the sample the low-mobility carriers cannot be moving in an 
aqueous olectrolytlc medium, which leaver only the possibility of such ions 
moving through the solid sample matrix (i.e. they travel through the minerals 
forming the rock). One may conclude, therefore, that the exponent ia I - I i ke decay 
(Figure 8a) is associated with properties of the minerals forming the rock. 

When exposing the sample to atmospheric moisture, after evacuation, the 
fast decay voltage Vp reappears in the rock response. This and related observations 
led the author (Alvarez, 1971) to suggest that, water molecules adsorbed on re-- 
sistive samples induce radical changes in their surface conductivity by forcing 
an electronic- type conduction in the sample surface. Consequently, in the model 
presently discussed, the fast decay voltage attributed to electronic conduction 
should bo related 1o surface conduction in the sample. 

One would expect fast relaxation times associated with electroni c-type con- 
duction, as well as a linear response of current vs voltage. The former feature 
is readily observed, although the term "fast relaxation" is taken here as relative 
to the relaxation time of the voltage V^, which is often of less than 10 seconds. 

The linear response seems to be supported by the linearity, above 3 volts, of 

the plot in Figure 8b, which coincides with the linear region of the l-V 
characteristic (Figure 8d). 

However, it would be more proper to decide on the linearity of Rp 
by obtaining the current associated only with Vp, which may be labeled Ip, and 
plot the lp"Vp characteristic. At the time the measurements were made we did 
not have means of obtaining Ip and thus the total current f-j. had to be used. 

The main differences between the Ip- Vp and the lj~ Vp plot* are expected to take 



place at voltages of values around the magnitude of the voltage shift from 
the origin Ci.e. around tl.Q yolt for the data in Figure 8d). 

Summarizing the discussion on the hematitic sandstone sample, and 
the schematic model of conduction proposed to explain its response, we have 
that: Cl) there is a linear region above 3 volts and a non-linear region 

below such a voltage; (2) a fast decaying voltage (Vp) is associated to 
surface conduction; (3) a slowly decaying voltage (V D ) is associated with 
conduction through the rock matrix, and (4) there is a persistent polariza- 
tion voltage obtained when the sample has not been excited, that corresponds 
to the displacement of the l-V character i stic from the origin. 

In addition, evidence of the presence of the polarization voltage 
mentioned in point (4) above was obtained when the sample was in a high 
vacuum: the currents flowing with forward and reversed polarities were 

different. From this fact follows that, whatever the mechanism responsible 
for such a polarization, it does not depend on the presence of water in the 
sample. 

The discussion of the resistive group of samples brings about a common 
feature present in all of them, namely the displacement of the l-V character- 
istics from the origin. We found that this feature may give rise to apparent 
non-linear behavior if R-V or £ plots of the data are made, but that 
the problem does not appear if the l-V characteristic is consider**!; from it, 
actual linear or non-linear behavior can be most properly inferred. 

The possibility of observing non-linear behavior at low voltages in 
rocks in contact with metallic electrodes is confirmed by the l-V plots of 
limestone, granite and hematitic sandstone. 

in the introductory discussion it was mentioned that non-linear behavior 
at low-voltages has been observed in relation to metal- insulator- metal or 



roetal-seraiconductor-roeta I systems, and tfvat space charge accumulations play 
an important role in the contact properties. 

By analogy with' such o&seryations we propose that, in resistivo rock 
samples In contact with metallic electrodes, the displacemenl of the l-V 
characteristics from the origin is due to space charge accumulations within 
the electrode-sample-electrode system. 

The spatial charge accumulations would give rise to the non-linearities 
in the case of the limestone, granite and herruititic sandstone samples. The 
linear behavior manifested by the argillaceous sandstone and the muscovite 
schist samples is not necessarily in disagreement with the presence of spatial 
charges since, as pointed out previously, contact effects can give rise to 
linear or non-linear behavior depending on the type of contact (ohmic, 
blocking or neutral) established between the metallic electrode and the 
sample, as well as on the magnitude of the applied voltage. Changes in the 
type of contact induced by the latter may be responsible for the transitions 
from non-linear to linear behavior in the granite and hematitic sandstone 
samples. 

Assuming that the displacement of the l-V characteristics from the 
origin is due to space charge accumulations, one can infer the existence of 

O 

the latter at high vacuum condi ti ns. from the evidence of polarization voltages 
obtained for the hematitic sandstone sample. This fact is of importance since 
it rules out the possibility of non-linearities arising from aqueous electro- 
lyte-electrode effects (i.e. the sample is extremely dry at high vacuum condi- 
tions) thereby supporting the hypothesis of the non-linearities being due 
to electrode-sample matrix effects. 



At any rate, electrode-sample effects are expected to be proportional to 
the magn Ltude of the difference between the sample conductivity and the elec- 
trodes' conductivity.^ From the data for galena and pyrite, presented below, 
one observes that the contrast in conductivity between the metallic electrodes 
and the sample is at least five orders of magnitude smaller for galena and 
pyrite than for the samples of the resistive group. Consequently, much smaller 
effects arising from the metal -sample interface are expected for the conductive 
samples than for the resistive ones. 

Figures 9a and b correspond to samples of the semiconductors galena and 
pyrite (i.e. "good"conductors) . As expected, the responses were linear. The 
experiments were not carried to higher voltages due to the low resistance of 
the samples which demanded larger than available currents; in any event the 
region of interest was the one around zero volts. Their responses to dc excita- 
tion voltages did not show the s low-decay region characteristic of resistive 
samples (Figure 8a). The fast decay Vp was equal to the applied excitation V q 
and therefore was zero for these samples; consequently the response to square 
pulses were square pulses for voltage and current. No variations of these 
responses were observed when the pyrite sample was evacuated to 10 ^ Torr. 

A formal proof of the hypothesis established above would require: (I) 
establishing the nature of the contact for each sample at different voltages 
and (.2) obtaining positive evidence of space charges in the electrode-sample 
interface. As of now we have been unable to determine the former; the latter 
is discussed in the next section. 

3 

The problem of the contrast in electrical properties between two media in 
relation to conduction phenomena has been previously analyzed (Alvarez, 1972). 
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SURFACE VOLTAGE DISTRIBUTION 


, The purposes of the_ experiment to be described next was to determine the 

effect of the metallic contacts on a resistive sample. If the properties 
of the metal-sample interface differ from those of the bulk of the sample, the 
voltage distribution in the electrode vicinity should manifest the difference. 

Figure 10 shows the schematic arrangement used in the experiment. The 
rock sample was a hematitic sandstone as the one used for the measurement 
reported in Figure 8. I ndi um-Mercury amalgam electrodes were used; one of them 
was constantly grounded while the other was polarized to + 4.9 volts; at times 
this electrode was floating or short-ci rcuited to ground depending on the 
requirements of our observations. A micro-electrode with a platinum tip exposed 
five microns was used to determine the voltage with respect to ground at 
points (X, Y) on the sample surface. All measurements were performed at dc 
•when transients had disappeared. Readings were taken at separations of tenths 
of a millimeter when close to the electrodes and by a millimeter elsewhere. Y was 
kept constant while X-travtrses were made between electrodes, then Y was 
Incremented by a millimeter and the process repeated. 

“The sample average dimensions were X = 24 mm, Y = 19 mm, and Z = 14 mm; 
the measurements were performed on 8 lines (i.e. Y = 0 through 7) covering 
somewhat less than half the surface of the sample. Figures I la, b, c, and d 
show the voltage values along four different lines; Y = 0 corresponds to the 
edge of the sample. 

The effect of the electrode-sample interface is manifested, at the positive 
electrode, by a voltage drop of around 1.7 volts occurring within a tenth of 

■v - * I 

a mi 1 1 imeter from the_eJectrode for Y = 0 and Y_= I.Ojron, ..and a voltage drop 



of 0.3 volts and 0.7 volts within the same distance occurs for Y = 4.0 mm and 
Y - 6.0 mm respectively. At the grounded electrode a similar effect is observed; 
'When going from the grounded electrode into the sample a "peak" is formed by 
a sharp increase followed by a decrease in voltage. Subsequently the average 
voltage rises when X decreases. 

From the data in Figures II a, b, c, and d, and recalling that those 
voltages correspond to stationary flow of direct current (i.e. the displacement 
current was zero when the measurements were taken), one can use the relation 
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to obtain Information about the charge distribution in the sample. The 
concentrations Cor charge densities) of positive and negative charge carriers 
at a given point are designated n + and n_; £ is the dielectric permittivity 

of the sample at the same point and e is the magnitude of the electron charge. 
Although there is insufficient information to apply the above relation 
properly (i.e. we do not know £ local ly ) we ^ ^ can draw some general 
conclusions from it. * 

The large voltage drop close to the positive electrode corresponds to an 
accumulation of negative charge in Its vicinity; that is n_ y n + in this 
region, while close to the grounded electrode a concentration of positive 
charge (I.e. n y n_) is observed. Along any line there is evidence of 
charge concentrations given by the increases and decreases in voltage over the 
average voltage slope. These reflect the inhomogeneous character of the 
sample, for in order to have steady state conditions in an inhomogeneous sample 
ft is necessary that charge accumulations occur at the interfaces between 
regions of different electrical properties (Alvarez, 1972). 



Since 
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T = _ (J-^V = coost. 

°x 2>x 


there Is obviously an Inverse relationship between the slope of the voltage 
and the value of the conductivity within two points separated by a distance 
A x. Therefore, associated to the space charge concentrations there are 
decreases in conductivity. The total current in the steady state was I = 1.9 x 

_5 

10 amp; the intensity of the electric field varies roughly from I volt/cm 

<b- 

to 170 volts/cm, thus variations of two orders of magnitude are observed 
between the average value of conductivity in the bulk of the sample and 
the regions close to the electrodes. 

In Figure 12, contours of the voltage values are presented; the lines 
show the equipotentials on the surface of the sample. The non-uniformity 
of the electric field Is evident. A word of caution is due regarding the 
contouring of the data: along Y-lines there are small vertical bars which 

show the actual points at which the measurements were taken; in order to 
•obtain the equipotentials a I i near variation of voltage was assumed between 
any two neighboring points. It is apparent that this is not a rigorous way 
of proceeding: suppose that, in the present case, measurements were taken 

—every other point and, assuming there is a linear variation of voltage 
between them, predict the values at intermediate points; clearly there will 
be places in which the predicted value and the actually determined one will 
show large differences. In spite of this Inaccurate assumption the general 
trends of voltage variation will manifest themselves when sufficient points are 
determined. In the case presently dealt with, closer sampling of the surface 
would modify the details of the contouring but the general trends would 


remain the same. 
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Earlier It was mentioned that deriving resistivity values from the total 
current, total voltage, and geometry of the sample in the case of measurements 
made with a two-electrode system (with parallel electrodes) was valid if 
the field E in the sample was uniform. The evidence presented above shows 
that assuming uniformity of the field may lead to a gross simplification of 
the actual conditions. Furthermore, it appears that by so doing essential 
information, that may help in explaining phenomena as the non-l Inearities, Is 
lost. 

Incidentally, the type Of resistivity values that would be obtained with 
a four-electrode system, having the same current electrodes and a separation 
of a few millimeters in the voltage electrodes, can be inferred from the 

data in Figure 12. The values of resistivity thus obtained evidently depend 
on position and would show scatter if various electrode separations were 
tried along the same line. 



CONCLUSIONS 


Non-I Inear behavior at voltages of less than 10 volts was observed for 
limestone, granite and hematitic sandstone samples in contact with metallic 
electrodes; argillaceous sandstone, muscovite schist, galena and pyrite 
showed linear behavior. AM of the resistive samples showed a displacement 
of the I —V characteristics from the origin, while the two samples of the 

jhwp 

con duct i ve^ di d not. 

Experimental determination of the voltage distribution on CL sample 
surface showed the existence of space charge concentrations on the electrodes’ 
vicinity. This fact strengthens further the hypothesis of non-linear behavior 
being due to contact effects. 

u. 

Evidence of space charge accumlation in a hematitic sandstone sample at 

A 

high vacuum conditions indicates that aqueous electrolyte-electrode effects 
are not the only possible cause of non-linearities. It was pointed out that 
interfacial phenomena were expected to be proportional to the magnitude of 
the difference between the sample conductivity and the electrodes’ conductivity. 

The use of dc excitation was suggested as the most effective way of 
observing electrode effects. Results obtained with a two-electrode system and 
dc excitation are relevant, in the context of our discussion, to measurements 
performed with four-electrode systems and ac. excitations. A brief discussion 
in this respect closes the section. 

in a four-electrode system there are contact effects between the sample 
and the electrolyte due to their different electrical properties ( i.e. the 
electrolyte js, the electrode). On the basis of the high conductivity contrast 



127 


between a sample and the electrolyte one may expect interfacial phenomena 
similar to the one reported in this work, in lieu of or in addition to 

i 

eletrochemica I phenomena occurring at the interfaces. 

Studies on synthetic ores using ac excitations have shown decreasing 
resistivities when the current density increases (McEuen et al, 1959; 

Scott and West, 1969), indicating the existence of non-linearities in rock- 
electrolyte systems. Furthermore, the presence of a bias voltage has been 
commonly observed by various researchers when measuring with ac excitations; 
it corresponds to the displacement of the l-V characteristics, at dc, that 
we have reported. In particular Scott and West, 1969, indicate that when 
measuring at (ac) equilibrium conditions a bias remained in their potential 
measuring circuit. They offer no explanation for the presence of such a bias, 
but it has to be attributed to electrolyte-sample effects since measurements 
with a four-electrode system preclude the possibility of electrode-sample 
effects. 

Thus, interfacial phenomena, non-linear effects and shifting of the l-V 
characteristics seem to be common features In twcKand four-^e lectrode systems, 
when either dc or ac excitations are used. 
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TABLE OF FIGURES 

Fig. I Linear and non-linear relations between (a) current vs voltage and 
(b) resistance vs voltage in a sample. 

Fig. 2 (a) R-V plot for a limestone sample. A non-linear region for V < 5 voth 

is obtained, while a linear region seems to start above 5 volts, (b) 

I - V characteristic of a similar sample of limestone; in this case one 
electrode was constantly grounded while the other underwent positive 
and negative variations in voltage. The non-linear region extends 
throughout the voltage region analyzed. Notice the displacement of 
the I —V characteristic from the origin, (c) The data in the previous 
figure is plotted as resistivity vs current density modulus. Since 
there are two different values of voltage for each current value in 
the I — V characteristic, there are two resistivity values for each value 
of the current density modulus. Zero resistivity corresponds to the 
point at which the characteristic interests the current axis, i.e. to the 
point at which there is non-zero current at zero voltage. 

Fig. 3. (a) A material assumed to be linear but having a displacement of its 

I — V characteristic from the origin is used to model the behavior of the 
limestone sample. The current and voltage scales are multiplied by 
arbitrary factors, (b) Assuming that the above material is made in a 
unit cube one can obtain the resistivity vs current density modulus plot. 
Non-linear behavior is shown for current density values smaller than 
a/!00J o , while for larger values a linear region is obtained. J 0 is 

the current flowing in the sample when its terminals are short circuited. 

+ 

The two resistivity branches ^ and y arise from the fact that the l-V 
characteristic does not go through the origin. 
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Fig. 4. 


Fig. 5. 
Fig. 6. 

Fig. 7. 


Fig. 8. 


(a) l-V characteristic of a granite sample. The non- I i near region 
appears below 3 volts, above this voltage the behavior is linear. The 
displacement from the origin is not appreciable in the scales used. 

(b) The correspond! ng R-V plot agrees with the l-V characteristic 

showing linear and non-linear behavior In the same voltage regions. 

The l-V plot of an argillaceous sandstone sample shows linear behavior; 
in the scales shown it appears not to be displaced from the origin. 

The l-V plot of a muscovite schist sample shows linear behavior, a 

displacement from the origin towards the positive side of the current - 

axis can be appreciated in the scales used. 

(a) R-V plot of the data in Figure 5 showing the effect of the dis- 
placement of the l-V characteristic towards the positive side of the 
current axis. It is manifested as a decreasing resistance when the 
zero voltage is approached from positive values. The points shown as 
a cross in a circle were obtained Immediately after the 9 volts 
excitation was applied to the sample. They indicate that an additional 
displacement of the l-V character! stic took place after application of 
the maximum voltage, (b) R-V plot of the data in Figure 6 gives no 
clear indication of linear behavior contrasting with the information 
given by the l-V characteristic. The point shown as a cross in a circle 
was obtained immediately after exciting the sample with 8 volts. However 
the sample recovers its original values if no new excitations are applied 
in several hours. 

(a) Response of a res i st i ve- type rock sample upon interruption of an 
externally applied voltage V c ; Vp and V^ are the magnitudes of the 
"fast" and "slow" decays, (b) Ratio of Vp to the total current in 
the sample, (c) Ratio of the total voltage to the total current in the 
sample, (d) Current vs voltage characteristic showing a displacement 
of + 1.0 volt from the origin. Non-linear behavior is shown below 3 
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Fig. 9. 
Fig. 10. 


Fig. II. 
Fig. 12. 


volts, above this voltage the characteristic is linear. 

R-V responses of (a) galena and (b) pyrite. Both show linear 
behavior. 

Schematic set-up for the determination of the voltage distribution 
in the surface of a sample. The micro-e lectrode traverses the 
lines taking voltage readings at points (X, Y) when a constant 
current is flowing. 

Voltage values as a function of distance X for a toematitic sandstone 
sample, (a) Y = 0, (b) Y = 1.0 mm, (c)*Y = 4.0 mm, and (d) Y = 6.0 mm. 
Equi potenti a Is on the surface of a hematitic sandstone sample showing 
the non-uniformity of the electric field when a constant current flows. 
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